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Several extensions of the Standard Model and in particular superstring theories 
suggest the existence of axion-like particles (ALPs), which are very light spin-zero 
bosons with a two-photon coupling. As a consequence, photon- ALP oscillations oc- 
cur in the presence of an external magnetic field, and ALPs can lead to observable 
effects on the measured photon spectrum of astrophysical sources. An intriguing 
situation arises when blazars are observed in the very-high-energy (VHE) band - 
namely above 100 GeV - as it is the case with the presently operating Imaging At- 
mospheric Cherenkov Telescopes (lACTs) H.E.S.S., MAGIC, CANGAROO III and 
VERITAS. The extragalactic background light (EBL) produced by galaxies during 
cosmic evolution gives rise to a source dimming which becomes important in the 
VHE band and increases with energy, since hard photons from a blazar scatter off 
soft EBL photons thereby disappearing into e'^e~ pairs. This dimming can be con- 
siderably reduced by photon- ALP oscillations, and since they are energy-independent 
the resulting blazar spectra become harder than expected. We consider throughout 
a scenario first proposed by De Angelis, Roncadelli and Mansutti - to be referred to 
as DARMA for short - in which the above strategy is implemented with photon- ALP 
oscillations triggered by large-scale magnetic fields, and we systematically investigate 
its implications for VHE blazars. We find that for ALPs lighter than 5 • 10~ eV the 
photon survival probability is larger than predicted by conventional physics above a 
few hundred GeV. Specifically, a boost factor of 10 can easily occur for sources at 
large distance and large energy, e.g. at 8TeV for the blazar lES 0347-121 at red- 
shift z = 0.188. This is a clear-cut prediction which can be tested with the planned 
Cherenkov Telescope Array and the HAWC water Cherenkov 7-ray observatory, and 
possibly with the currently operating lACTs as well as with detectors like ARGO- 
YBJ and MILAGRO. Moreover, we show that the DARMA scenario offers a new 
interpretation of the VHE blazars detected so far, according to which the large spread 
in the values of the observed spectral index is mainly due to the wide spread in the 
source distances rather than to large variations of their internal physical properties. 
Finally, we stress that ALPs with the right properties to produce the above effects 
can be discovered by the GammeV experiment at FERMILAB and more likely with 
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the planned photon regeneration experiment ALPS at DESY. 

PACS numbers: 14.80.Mz, 95.30.-k, 95.85.Pw, 95.85.Ry, 98.70.Rz, 98.70.Vc, 98.70.Sa 



I. INTRODUCTION 

The Standard Model (SM) of strong, weak and electromagnetic interactions with three 
sequential families of quarks and leptons provides a satisfactory description of the experi- 
mental results concerning elementary-particle physics concerning elementary-particle physics 
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at energies up to about the Fermi scale of weak interactions Gp ~ 250 GeV. 

Yet, nobody would seriously regard the SM as the ultimate theory of fundamental pro- 
cesses. Apart from more or less aesthetic reasons and actual problems like the g — 2 value of 
the muon and the top forward-backward asymmetry such an expectation is made compelling 
by the observational evidence for non-baryonic dark matter ultimately responsible for the 
formation of structure in the Universe as well as for dark energy presumably triggering the 
present accelerated cosmic expansion. 

So, the SM is presently viewed as the low-energy manifestation of some more fundamental 
and complete theory of all elementary-particle interactions including gravity. Every specific 
approach to extend the SM in such a way is characterized by a set of new particles along 
with their mass spectrum and their interactions with the standard world. 

Although it is presently impossible to tell which proposal out of so many ones has any 
chance to successfully describe Nature, it looks remarkable that attempts along very different 
directions such as four-dimensional ordinary and supersymmetric models [Ij, Kaluza-Klein 
theories |21 13] and especially superstring theories [IHS] all suggest the existence of axion-like 
particles (ALPs) [7]. A general argument supporting this conclusion will be given in Section 
II. 

ALPs are very light pseudo-scalar spin-zero bosons a characterized by a two-photon 
coupling 077. As the name itself suggests, they are a sort of generalization of the axion, the 
pseudo-Goldstone boson associated with the Peccei-Quinn symmetry proposed as a natural 
solution to the strong CP problem PHTT]. But while the axion enjoys a strict relationship 
between its mass and the 077 coupling constant, these two parameters are to be regarded 
as unrelated for ALPs. In fact, depending on the actual values of their mass and 077 
coupling constant, ALPs can play an important role in cosmology, either as cold dark matter 
particles [I2] or as quintessential dark energy [T3] . 

A remarkable consequence of the 077 coupling is the phenomenon of photon- ALP mixing, 
which takes place in the presence of an external electromagnetic field and leads to two distinct 
effects. One is photon- ALP oscillations [131 [15] . which is quite similar to the oscillations of 
massive neutrinos with different flavours. The other consists in the change of the polarization 
state of photons traveling in a magnetic field [ini ITU] . 

It turns out that ALPs are extremely elusive in high-energy experiments and the only 
way to look for them in the laboratory requires very careful polarimetric measurements to 
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be carried out on a laser beam [T7] or alternatively a photon regeneration experiment to 
be performed [18]. Successful detection of ALPs in present-day experiments of this kind is 
possible in either case for a fairly large 077 coupling. 

Astrophysical manifestations of ALPs appear the best way to discover their existence, 
since they can give rise to observable effects even for values of the 077 coupling constant 
much smaller than those tested so far in laboratory experiments. Indeed, it is known since 
a long time that for values of the 077 coupling constant that look hopelessly small to be 
probed today in the laboratory the stellar evolution would be dramatically altered [19] and 
this fact sets a strong upper bound on the coupling in question, which is consistent with the 
negative result of the CAST experiment at CERN [20] . 

In the last few years it has been realized that photon-ALP oscillations triggered by 
intervening cosmic magnetic fields along the line of sight can produce detectable effects in 
observations of bright X-ray and 7-ray sources [^TH^ . The effect becomes larger as the 
distance of the sources increases: blazars, which are Active Galactic Nuclei (AGN) with 
the beam pointing towards us, constitute the most distant long-lasting gamma-ray sources 
observed, and are thus the most obvious case study. In order to bring out most simply the 
relevance of ALPs in the present context, we neglect cosmological effects at this introductory 
level. 

The mean free path of very-high-energy (VHE) photons - namely with energy above 
100 GeV - is limited by their interaction with background photons in the Universe through 
the process 77 — )• e^e~. A high energy photon from a distant blazar has a nonnegligible 
probability to scatter off background photons in the optical/near infrared band permeating 
the Universe - the so-called extragalactic background light (EBL) - thereby disappearing into 
an e~^e~ pair [321 ES]- The VHE photon mean free path depends on the EBL density. Re- 
cently it has become possible to model accurately such density: several parametrizations are 
available, essentially coincident, and throughout this paper we employ the recent EBL model 
of Franceschini, Rodighiero and Vaccari (FRV) [M])- Figure [l] shows the pair-production 
mean free path of a VHE photon as a function of its energy E within the FRV model. 

The effects of the photon-ALP oscillations on the gamma yield from blazars can be 
summarized as follows: 

• For E < 100 GeV we infer from Figure [l] that A-y(i?) is comparable with the Hubble 
radius, and so EBL absorption is negligible. In such a situation photon-ALP oscil- 
lations can only give rise to a source dimming above a certain energy threshold E^,. 
Hence a characteristic distortion of the source spectrum around E^, is the observable 
prediction, which can be searched for with the Fermi/LAT mission if i?^, happens to 
lie in the instrument energy range, namely for 30MeV < E^ < 300 GeV [23] . 

• For E > 100 GeV Figure [l] shows that EBL effects become important since now A^(£') 
quickly decreases as E increases. Once emitted, photons can convert into ALPs and 
next reconvert back into photons before reaching the Earth. A possibility is that 
photon-ALP oscillations take place in intergalactic space and the resulting scenario 
has been called DARMA (acronym for De Angelis, Roncadelli and Mansutti) [231 - I26] . 
Alternatively, the 7 — )• a conversion can occur inside the blazar while the a —)■ 7 re- 
conversion can happen in the Milky Way [2?]. Of course, also both options can be 
realized [30]. In all cases, photons acquire a split identity, travelling for some time as 
real photons and for some time as ALPs. However, they suffer EBL absorption only 
when they are real photons, which means that the effective photon mean free path 
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FIG. 1: The pair-production mean free path of a VHE photon is plotted versus its energy 
E within the EBL model of FRV. Only conventional physics is assumed and in particular the 
possibility of photon- ALP oscillations is ignored. 

A^^eff(-S) is actually larger than X^{E) as predicted by conventional physics (see Figure 
[T]). Since the photon survival probability depends exponentially on minus the optical 
depth - which in turn goes like the source distance divided by Xy^es{E) - even a slight 
increase of A^^eff(-E) with respect to X-y{E) produces a substantial enhancement of the 
photon survival probability and so of the observed flux. Note that for a given detector 
sensitivity a larger photon survival probability just means that a larger distance can be 
probed, so that the VHE Universe becomes more transparent than generally believed. 
Actually, since the EBL absorption increases with energy whereas the photon- ALP os- 
cillation probability is energy- independent, the observed flux enhancement gets larger 
and larger as the energy increases. As a consequence, the observed spectra are harder 
than currently expected. 

Thus - depending on the values of the free parameters - a hardening of the observed blazar 
spectra is the main prediction of photon-ALP oscillations concerning the VHE band between 
100 GeV and 100 TeV, which can be probed by Imaging Atmospheric Cherenkov Telescopes 
(lACTs). More in detail, the presently operating lACTs H.E.S.S., MAGIC, CANGAROO 
III and VERITAS can reach up ~ 20 TeV with difficulty, whereas the planned Cherenkov 
Telescope Array (CTA) and the HAWC water Cherenkov 7-ray observatory will be able to 
explore the whole VHE band with a much larger sensitivity. We should add that the VHE 
range can also be analyzed by other available detectors, like the Extensive Air Shower arrays 
ARGO-YBJ and MILAGRO. 

The aim of the present paper is to investigate the DARMA scenario in great detail by 
employing the EBL model of FRV [31j, and systematically working out all its implications 
for VHE blazar observations. Everything is calculated starting from first principles. 

More specifically, our primary goal is to evaluate the photon survival probability within 
the DARMA scenario, which allows in turn to quantify the resulting hardening of VHE 
blazar spectra and to determine how much the Universe becomes more transparent to VHE 
photons than usually thought. 
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We find that these effects can be quite substantial for the free parameters in allowed 
realistic ranges provided that the ALP mass satisfies the condition m < 5-10~^° eV. Actually, 
our prediction can become spectacular above an energy threshold which is well below the 
upper detection limit of the CTA and the HAWC observatory. Denoting by £"10 the energy at 
which the photon survival probability is 10 times larger than that dictated by conventional 
physics, £"10 turns out to decrease with the source redshift z. In the most favourable case, 
we get Eio - 30 TeV for Mrk 421 at z = 0.030, Eio ~ 8 TeV for lES 0347-121 at z = 0.188, 
Eio ^ 2 TeV for 3C 66A at ;z = 0.444 and Eio ^ 2 TeV for 3C 279 at z = 0.536. In addition, 
it follows from the above considerations that for energies larger than Eio the photon survival 
probability gets even more enhanced. 

In spite of the fact that only future observations can provide a clear-cut check of the 
DARMA scenario - and so ultimately of the existence of an ALP lighter than 5 ■ 10~^° eV - 
it looks natural to inquire whether available data from lACTs contain some hints in favour 
of this scenario. As we shall see, the emitted as well as the observed spectra of VHE blazars 
detected so far have to a good approximation a single power-law behaviour [35] , so that they 
are characterized by their slope, which - up to a minus sign - is the emitted Fem and observed 
Fobs spectral index, respectively. So, we can rephrase the above conclusion by stating that 
for a fixed value of Fgm pertaining to a given source the DARMA scenario predicts that 
under suitable conditions Fobs should be smaller than within conventional physics. 

In 2006 the H.E.S.S. collaboration reported some evidence that the spectra of the two 
blazars H 2356-309 and lES 1101-232 have Fobs smaller than expected and this fact was 
interpreted as strongly suggesting an EBL attenuation lower than currently believed [35]. A 
similar conclusion emerged in 2007 with the discovery of the blazar 3C 279 by the MAGIC 
collaboration [37] . Based on preliminary data [3H] and a specific EBL model [39j, in a 
previous paper [21] two of us (A. D. A. and M. R. together with O. Mansutti) have shown 
for the first time that the mechanism of photon- ALP oscillations can substantially reduce 
the EBL attenuation for distant blazars and in particular that it can successfully explain 
the observed spectrum of 3C 279 for allowed realistic values of the free parameters. 

Subsequent developments have demonstrated that realistic EBL models account for VHE 
blazar observations without the need of any unconventional physics, provided that the large 
spread in the values of Fobs is fully traced back to an equally large spread in the values of 
Fem- Further, far-away sources - for which EBL absorption is a large effect - turn out to 
have energy spectra similar to those of some nearby blazars, for which EBL attenuation is 
negligible. This means that for distant sources Fem has to be considerably smaller than for 
nearby ones. Even though a physical explanation for the occurrence of very small values of 
Fem has recently been achieved [101 - H2] . one is nevertheless led to the cosmic opacity problem, 
namely to wonder why these physical effects are important for distant blazars only. 

We show that within the DARMA scenario the situation is quite different. As a result 
of the competition between EBL attenuation and photon-ALP oscillations, two important 
conclusions emerge: 

• The values of Fgm for far-away VHE blazars are in the same ballpark of nearby ones, 
so that the cosmic opacity problem is solved. 

• The observed large spread in the values of Fobs arises mainly from the wide spread in 
the source distances while the required scatter in the values of Fgm is small. 

Finally, we stress that ALPs with the right properties to produce the above effects can be 
discovered by the GammeV [13] experiment at FERMILAB and more likely by the planned 
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photon regeneration experiment ALPS at DESY [H] or with large xenon scintillation detec- 
tors developed for dark matter searches |15] . Thus, it looks amazing that the discovery of a 
new particle - besides very important in its own right - would also provide a sort of glasses 
that allow us to watch much farther out into the 7-ray Universe. 

The paper is organized as follows. In Section II we review the motivation and the main 
properties of ALPs. Particular attention is paid to the propagation of a photon/ALP beam 
in the presence of a homogeneous magnetic field, a problem that can be solved exactly. 
Also the astrophysical and cosmological constraints on the ALP parameters are considered. 
Section III addresses in great detail the properties of VHE blazar spectra. Everything in this 
Section is discussed within conventional physics and at this stage ALPs are totally neglected. 
After a brief account of cosmic opacity from a microscopic point of view, various models of 
the EBL are cursorily described and compared. Then it is shown that conventional physics 
provides a correct qualitative understanding of the observed blazar spectra but at the same 
time it requires some degree of correlation between Pem and z which gives rise to the cosmic 
opacity problem. It is also argued that a way out of this problem calls for an unconventional 
photon propagation, and various proposals are briefly reviewed. In Section IV the main 
conclusions drawn in Sections II and III are combined together to build up the DARMA 
scenario, which is investigated in detail. The phenomenon of photon-ALP oscillations is 
considered here in the presence of EBL absorption as well as for a domain-like configuration 
typical of large-scale magnetic fields. The photon survival probability is ultimately computed 
by a numerical code and it arises as an average over 5000 random realizations of the beam 
propagation from the source to us, each corresponding to randomly chosen directions of the 
magnetic field inside every domain. Sections V, VI and VII are devoted to the discussion of 
the implications of the DARMA scenario for VHE blazar observations. A detailed analysis 
of the behaviour of the photon survival probability for some representative values of the free 
parameters is presented in Section VI, where its relevance for future VHE blazar observations 
is stressed. Section VII offers a new interpretation of the observed VHE blazars, which solves 
the cosmic opacity problem and traces the large spread in the values of Pobs niainly to the 
wide spread in the source distances while the demanded scatter in the values of Pem is small. 
Finally, we offer our conclusions in Section VIII. A convenient method to solve the eigenvalue 
problem for a 2 by 2 matrix with complex coefficients is presented in Appendix A, whereas 
an approximate analytic evaluation of the optical depth within the FRV model of the EBL 
is reported in Appendix B. 

II. AXION-LIKE PARTICLES (ALPS) 

We review the conceptual motivations in favour of ALPs as well as their properties that 
are most relevant for our further needs. Natural Lorentz-Heaviside units with h = c = kB = 1 
are employed throughout the paper unless otherwise stated. 

A. Motivation 

As already stressed, the Standard Model (SM) of particle physics is presently regarded 
as the low-energy manifestation of some more fundamental theory (FT) characterized by a 
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very large energy scale A ^ Gp , with Gp ~ 250 GeV. We collectively denote by </> 
the SM particles together with possibly new undetected particles with mass smaller than 
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Gp^^'^, while all particles much heavier than Gp^^'^ that are present in the FT are collectively 
represented by $. Correspondingly, the FT is defined by a Lagrangian of the form £ft(0, ^) 
and the generating functional for the corresponding Green's functions reads 

Zft[J,K] = N j j V<!> exp (^i j d^x £ft(0,$) + (^J + <^K ^ , (1) 

where J and K are external sources and is a normalization constant. The resulting low- 
energy effective theory then emerges by integrating out the heavy particles in ZptIJjK], 
and so the low-energy effective Lagrangian Ces{4') is defined by 

exp(^i j d^xC^s{(f))^ = j V<i>exp(^i j rf^x£FT(0,$)^ • (2) 

Evidently, the SM Lagrangian is contained in £efj(0), and - in the absence of any new 
physics below Gp - it will differ from Ces{(p) only by non-renormalizable terms involving 
the (p particles alone, that are suppressed by inverse powers of A. 

In any theory with a sufficiently rich gauge structure - which is certainly the case of 
the FT - some global symmetry Q invariably shows up as an accidental consequence of 
gauge invariance. Since the Higgs fields which spontaneously break gauge symmetries carry 
nontrivial global quantum numbers, it follows that the group Q undergoes spontaneous sym- 
metry breaking as well. As a consequence, some Goldstone bosons - which are collectively 
denoted by a if ^ is non-abelian - are expected to appear in the physical spectrum and their 
interactions are described by the low-energy effective Lagrangian, in spite of the fact that 
Q is an invariance of the FT. We stress that Goldstone bosons are necessarily pseudo-scalar 
particles |16]. 

As far as our main line of development is concerned, the FT is supposed to describe 
quantum gravitational effects and it is a common lore that they always explicitly break 
global symmetries In fact, this point can be understood in an intuitive fashion. Since 
black holes do not possess any definite global charges, global symmetries are violated in 
any scattering process involving black holes. So, we end up with the general conclusion 
that provided that the Lagrangian of the FT possesses some spontaneously broken global 
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symmetry then pseudo- Goldstone bosons with mass much smaller than Gp are necessarily 
present in the low-energy effective Lagrangian. 

Therefore, by splitting up the set into the set of SM particles 0sm plus the pseudo- 
Goldstone bosons a, the low-energy effective Lagrangian has the structure 

^es{(pSM, 0.) = -^^SmI^Sm) + 'Cnonren(0SM) + Acn(a) + Aen(0SM, + -^nonrenl^SM, 0.) , (3) 

where jCren{4>SM, d) stands for renormalizable soft-breaking terms that can be present when- 
ever Q is not an automatic symmetry of the low-energy effective theory [48j . 

— 1/2 

Needless to say, it can well happen that between Gp and A other relevant mass scales 
Ai, A2, ... exists. In such a situation the above scheme remains true, but then Q may be 
spontaneously broken at such an intermediate scale. 

Finally, we would like to stress that a very thoroughly analysis by Arvanitaki et al. [6] in 
the context of superstring theories and by Turok [3] in fundamental theories with compact 
extra dimensions have made the above conclusion more specific, showing that in either case 
the pseudo-Goldstone bosons are actually ALPs. 
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B. Axion as a prototype 

A characteristic feature of the SM is that non-perturbative effects produce the term 
A£e = 6'(7^G'^''Ga^j,/327r^ in the QCD Lagrangian, where 9 is an angle, g and are 
the gauge couphng constant and the gauge field strength of S'f/c(3), respectively, and 
G^" = \e^^'"^Gapa- AH values of 6 are allowed and theoretically on the same footing, but 
nonvanishing 9 values produce a P and CP violation in the strong sector of the SM. An 
additional source of CP violation comes from the chiral transformation needed to bring 
the quark mass matrix M.q into diagonal form, and so the total strong CP violation is 
parametrized by ^ = + arg Det A^^. Observationally, a nonvanishing 9 would show up 
in an electric dipole moment dn for the neutron. Consistency with the experimental upper 
bound \dn\ < 3 ■ 10~^^ e cm requires |^| < 10~^ Thus, the question arises as to why \9\ is 
so unexpectedly small. A natural way out of this fine-tuning problem - which is the strong 
CP problem - was proposed by Peccei and Quinn [S] over 30 years ago. Basically, the idea 
is to make the SM Lagrangian invariant under an additional global f/(l)pQ symmetry in 
such a way that the A^g term can be rotated away. While this strategy can be successfully 
implemented, it turns out that the [/(1)pq is spontaneously broken and then a Goldstone 
boson is necessarily present in the physical spectrum. Things are slightly more complicated, 
because [/(1)pq is also explicitly broken by the same non-perturbative effects which give rise 
to A£0. Therefore, the would-be Goldstone boson becomes a pseudo- Goldstone boson - the 
original axion [H] - with nonvanishing mass given by 

m~0.6f — J J eV , (4) 

where denotes the scale at which f/(l)pQ is spontaneously broken. Qualitatively, the 
axion is quite similar to the pion and it posseses Yukawa couplings to quarks which go like 
the inverse of fa- Moreover - just like for the pion - a two-photon coupling 077 of the 
axion a is generated at one- loop via the triangle graph with internal fermion lines, which is 
described by the effective Lagrangian 

'^a77 = -T^i^^'^^M^a = ^E-Ba , (5) 

0-77 iiu ^ ■, \ J 

where Ff^" = (E, B) = d^A'^ — d'^A^^ is the usual electromagnetic field strength and F^'^ = 
^e^'^P^'Fpf^. The constant M entering Eq. (p| has the dimension of an energy and is given by 



M=1.2-lO''k(^J^] GeV, (6) 



with k a model-dependent parameter of order one |19]. Note that M oc fa and turns out 
to be independent of the mass of the fermions running in the loop. Hence, the axion is 
characterized by a strict relation between its mass and two-photon coupling 

,„^0.7*(15^)eV. (7) 

In the original proposal [8], [/(1)pq is spontaneously broken by two Higgs doublets which 

— 1/2 

also break SU^a{2) x Uy{l) spontaneously, so that fa < Gp . Correspondingly, from eq. 
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Q we get m > 24KeV. In addition, the axion is rather strongly coupled to quarks and 
induces observable nuclear de-excitation effects [50]. In fact, it was soon realized that the 
original axion was experimentally ruled out ^T\. 

A slight change in perspective led shortly thereafter to the resurrection of the axion 
strategy. Conflict with experiment arises because the original axion is too strongly coupled 
and too massive. But, given the fact that both m and all axion couplings go like the inverse 
of fa the axion becomes weakly coupled and sufficiently light provided that one arranges 

— 1/2 

fa ^ Gp . This is straightforwardly achieved by performing the spontaneous breakdown 
of [/(l)pQ with a Higgs field which is a singlet under SUu,{2) x Uy{l) [lO]. Note that we are 
thereby led to the conclusion that the f/(l)pQ symmetry has nothing to do with the low- 
energy effective theory to which the axion belongs, but rather it arises within an underlying 
more fundamental theory. 

Thus, we see that the axion strategy provides a particular realization of the general 
scenario outlined in Subsection II-A, with Q = [/(1)pq, Ai = fa and i3nonren(0SM, including 
£a77 among other terms involving the SM fermions. This fact also entails that new physics 
should lurk around the scale at which ?7(1)pq is spontaneously broken. The same conclusion 
is reached from the recognition that the Peccei-Quinn symmetry is dramatically unstable 
against a tiny perturbation - even at the Planck scale - unless it is protected by some 
discrete gauge symmetry which can only arise in a more fundamental theory [52|. 



C. Beyond the axion: ALPs 

A generic feature of many extensions of the SM along the lines discussed in Section I is 
the prediction of ALPs. Generally speaking, ALPs are a straightforward generalization of 
the axion but important differences exist between the axion and ALPs mainly because the 
axion arises in a very specific context while in dealing with ALPs the aim is to bring out 
their properties in a model- independent fashion as much as possible [7]. This attitude has 
two main consequences: 

• Only ALP-photon interaction terms are taken into account. Therefore, any other 
possible coupling of ALPs to SM particles is presently discarded and this entails that 
^nonren(</'SM5 a) in Eq. ^ Only iucludcs Ca-y-y as defined by Eq. (|5]). Observe that 
such an ALP coupling to two photons 077 is just supposed to exist without further 
worrying about its origin. 

• The parameters m and M are to be regarded as unrelated for ALPs, and it is merely 

— 1 /2 ~ 1 /2 

assumed that m <^ Gp and M ^ Gp 
As a result, ALPs are described by the Lagrangian 

£alp = - d^a d.a --m^a^- FauF^^a = - 9„a - - + E ■ B a. (8) 

2 ^2 AM ^ 2^2 M ^ ' 

D. Photon- ALP mixing 

What ultimately characterizes ALPs is the trilinear 077 vertex in £alp, which gives rise 
to photon-ALP mixing in the presence of an external magnetic field B. More specifically, 
what happens can be described as follows. 
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In such a situation, an off-diagonal element in the mass matrix for the photon-ALP system 
shows up. Therefore, the interaction eigenstates differ from the propagation eigenstates and 
the phenomenon of photon-ALP oscillations shows up [T^ US] . This is analogous to what 
takes place in the case of massive neutrinos with different flavours, apart from an important 
difference. All neutrinos have equal spin, and so neutrino oscillations can freely occur. 
Instead, ALPs are supposed to have spin zero whereas the photon has spin one, hence one 
of them can transform into the other only if the spin mismatch is compensated for by an 
external magnetic field. Note that the strength of this effect depends on the ratio B/M and 
not on B and M separately. 

We denote by E the electric field and by k the wave vector of a propagating photon 
at a given space-time point. Further, let and be the components of the external 
magnetic field B along k and perpendicular to k, respectively. Because E is orthogonal 
to k it follows that only the term E ■ B-y survives in £alp- We next split up E into two 
components, one Ey in the plane defined by k and B and the other E^ perpendicular to 
that plane. By construction, E_l is orthogonal to B-r, and so the 077 coupling in £alp goes 
like E\\ Bxa, which exhibits two characteristic properties of ALPs. First, the photon-ALP 
mixing depends only on the transverse component By of the external magnetic field; for 
notational simplicity we will write B rather than Bt in the following. Second, only photons 
linearly polarized along Ey actually mix with ALPs, whereas photons with polarization Ej^ 
do not mix. As a consequence, the 077 coupling acts like a polarimeter, in the sense that it 
gives rise to a change of the photon polarization state. This effect can be used to look for 
ALPs both in high-precision polarimetric measurements performed in the laboratory [T7j 
and in certain astrophysical observations in which the polarization state of the detected 
photons can be measured [3T] . 



E. Photon/ ALP beam propagation 

We shall be concerned throughout with a monochromatic, unpolarized photon/ALP beam 
of energy E and wave vector k propagating in a cold medium which is both magnetized 
and ionized (from now on E denotes the energy, and since the electric field will never 
be considered again no confusion arises). We suppose for the moment that the external 
magnetic field B is homogeneous and we denote by rig the electron number density. We 
employ an orthogonal reference frame with the y-axis along k, while the x and z axes are 
chosen arbitrarily. 

It can be shown that in this case the beam propagation equation following from £alp 
can be written as [12] 

+ + 2^-^o) m = (9) 

with 

/ My) \ 

ij{y) ^ AM , (10) 

where Ax{y) and Az{y) denote the photon amplitudes with polarization (electric field) along 
the X- and 2;-axis, respectively, while a{y) is the amplitude associated with the ALP. It is 
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useful to introduce the basis {|7x), 17^), !«)} defined by 

l7.) = 1^ j , (11) 

|7.) = f 1 I , (12) 



\a)=\^O^j , (13) 

where I'jx) and |7^) represent the two photon hnear polarization states along the x- and 
2;-axis, respectively, and \a) denotes the ALP state. Accordingly, we can rewrite ip{y) as 

tP{y) = Ax{y) |7,) + A,{y) I7,) + a(y) \a) , (14) 

and the real, symmetric photon-ALP mixing matrix Aio entering Eq. ^ has the form 

^xx ^xz 

Mo = I A.. A.. I , (15) 

Al, Af„^ A. 



07 a7 



where we have set 



A"^ = — fl6) 
A^ = — fl7l 

A„. = - ^ . (18) 

While the terms appearing in the third row and column of Aio are dictated by £alp 
and have an evident physical meaning, the other A-terms require some explanation. They 
reflect the properties of the medium - which are not included in £alp - and the off-diagonal 
A-terms directly mix the photon polarization states giving rise to Faraday rotation. 

In the present paper we are interested in the situation where the photon/ ALP energy is 
much larger than the ALP mass, namely E ^ m. As a consequence, the short-wavelength 
approximation can be successfully employed and can be implemented as [15] 

f <P f d \ f d \ , s. f d 



which turns the second-order beam propagation equation ^ into the first-order one 

(^i^ + E + Mo^ i'{y)=0. (20) 



We see that a remarkable picture emerges, wherein the beam looks formally like a three- 
state nonrelativistic quantum system. Explicitly, they are the two photon polarization states 
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and the ALP state. The evolution of the pure beam states is then described by the three- 
dimensional wave function ip{y) - with the y-coordinate replacing time - which obeys the 
Schodinger-like equation ( [20| ) with Hamiltonian 

Ho = -{E + Mo) . (21) 



Denoting by Uo{y,yo) the transfer matrix - namely the solution of Eq. (20) with initial 
condition Uo{yo,yo) = 1 - the propagation of a generic wave function can be represented as 

^{y) = Uo{y,yo)^{yo) . (22) 

Moreover, we have 

Uo{y,yo) = e'''^y'y^^^Uo{y,yo) , (23) 
where Uo{y,yo) is the transfer matrix associated with the reduced Schodinger-like equation 

%^^') '^(^) = °- (2^) 

Because B is supposed to be homogeneous, we have the freedom to choose the z-axis along 
B, so that Bx = 0. The diagonal A-terms receive in principle two different contributions. 
One comes from QED vacuum polarization, but since we will be dealing with very weak 
magnetic fields this effect is irrelevant [15] . The other contribution arises from the fact that 
the beam is supposed to propagate in a cold plasma, where charge screening produces an 
effective photon mass resulting in the plasma frequency 

where a is the fine-structure constant and rrig denotes the electron mass, which entails 



Ap. = -^. (26) 

Finally, the A^.^, A^x terms account for Faraday rotation, but since we are going to take 
E in the VHE 7-ray band Faraday rotation is negligible. Altogether, the mixing matrix 
becomes 

/Api \ 
M^o^ = Api A,^ , (27) 

\ A„^ AaaJ 

with the superscript (0) recalling the present choice of the coordinate system and 

A„ - ^ . (28) 

We see that decouples away while only mixes with a, showing that in the present 
approximation plasma effects do not change the qualitative features previously found in 
vacuo. 

Application of the discussion reported in Appendix A with Ai — t- A^q^'' yields for the 
corresponding eigenvalues 

Ao,i = Api , (29) 
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A, 



0,2 



A, 



0,3 



- Api + A,, + A 

OSC / ) 



(30) 
(31) 



where we have set 



A. 



(Api - Aaa)' + 4 (A 



nl/2 



a7 ) 



m — u 



2 \ 2 

pi 



2E 



B 
M 



1/2 



(32) 



As a consequence, the transfer matrix associated with Eq. (24) with mixing matrix Ai 



can be written with the help of Eq. (165) as 



,(0) 



(33) 



where the matrices To^i(O), To^2(0) and To^3(0) are just those defined by Eqs. (166), ( |167 ) 
and (168) as speciahzed to the present situation. Actually, a simplification is brought about 



by introducing the photon- ALP mixing angle 



a = - arctg 
2 ^ 



2A 



ay 



arctg 



2E 



B 



M J \ m? — 



(34) 



since then simple trigonometric manipulations allow us to express the above matrices in the 
simpler form 

'10 

To,i(0) = I I , (35) 




^0,3(0) 





sin^ a — sin a cos a 

— sin a cos a cos^ a 



cos^ a sin a cos a 
sin a cos a sin^ a 



(36) 



(37) 



Now, the probability that a photon polarized along the 2;-axis oscillates into an ALP after 
a distance y is evidently 

^i°U(z/) = l(a|Wo(l/,0;0)|7.)|' (38) 
and in complete analogy with the case of neutrino oscillations [19] it reads 



PoX^aiy) = sin^2a sin^ ^ 



(39) 



which shows that Aqsc plays the role of oscillation wave number, thereby implying that the 



oscillation length is Lose = 27r/Aosc- Owing to Eq. (34), Eq. (39) can be rewritten as 

2 



(0) 



0,7z 



B 



MAr 



sm 



AoscZ/ 



(40) 
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which shows that the photon-ALP oscillation probability becomes both maximal and energy- 
independent for 

Aosc ^ ^ , (41) 



M 

and explicitly reads 

P^l^Ay) - sin^ ( ^ ) . (42) 



This is the strong-mixing regime^ which - from the comparison of Eqs. (32) and (41) - turns 
out to be characterized by the condition 



\m^-ujl^\ B 



« T7 > (43) 



2E M 

and so it sets in sufficiently above the energy threshold 

E. ^ . (44) 

Note that the strong-mixing regime is unbounded from above since the contribution to 
photon-ALP mixing arising from QED vacuum polarization is negligible for the magnetic 
fields considered in this paper |29j . 

Below the photon-ALP oscillation probability becomes energy-dependent and vanish- 
ingly small. 

So far, our discussion was confined to the case in which the beam is in a pure polarization 
state. This assumption possesses the advantage of making the resulting equations partic- 
ularly transparent but it has the drawback that it is too restrictive for our analysis. For, 
photon polarization cannot be measured in the VHE 7-ray band, and so we have to treat the 
beam as unpolarized. As a consequence, it will be described by a generalized polarization 
density matrix 

/ My) \ 

p{y) = My) ® ( My) My) a{y) )* (45) 

V a{y) J 

rather than by a wave function ipiy)- Remarkably, the analogy with non-relativistic quantum 
mechanics entails that p{y) obeys the Von Neumann-like equation 

z| = [p,A^o] (46) 



associated with Eq. (24). Thus, the propagation of a generic p{y) is given by 

p{y) =Uo{y,yo) p{yo)l^l{y,yo) (47) 

and the probability that a photon/ ALP beam initially in the state pi will be found in the 
state p2 after a distance y is 

Po,p,^pM = Tr(p2Wo(z/,0)piWo^(i/,0)) , (48) 

since we are assuming as usual that Trpi = Trp2 = 1. Observe that in Eqs. (46), (47) 
and (48) we have dropped the superscript (0) in Aio and replaced Uo{y,yo]0) by Uo{y,yo) 
because they retain their form for an arbitrary choice of the coordinate system. 
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In view of our subsequent discussion it proves essential to deal with the general case in 
which B is not aligned with the z-axis but forms a nonvanishing angle ip with it. Corre- 
spondingly, the mixing matrix Ai^ presently arises from Ai^^^ through the similarity trans- 
formation 

Mo = V\ij)M''oW{ij) (49) 
operated by the rotation matrix in the x-z plane, namely 



cos ip — sin ip 
y (■?/') =1 smijj cosip 
1 



(50) 



This leads to |28] 



Mr 









Api 



Aa'y sinip Aa-y cosip 



Aa-y COSIp 



(51) 



indeed in agreement with Eq. (15) within the considered approximation. Therefore the 
transfer matrix reads 

Uo{y,yo;i^) = V\ij)Uo{y,yo;0)V{ij) 
and its explicit representation turns out to be 

l^oiy, yo; V^) = e'''^y~y'^ To,i(7/>) + e^^^(^-^») To,2W + e^^«(^-^°) To,3(^) , 



with 



cos^ -ip — sin tjj cos -ip 
To^i{ip) = j — sin -i/; cos ^/^ siii^ ip 




To,3(^) 



sin^ 9 sin^ ip sin^ a sin ip cos ip — sin a cos a sin ip 
sin^ a sin ip cos ip sin^ a cos^ ip — sin a cos a cos ip 
— sin a cos a sin ^/^ — sin a cos a cos ^/^ 



cos^ a 



sin^ i/' cos^ a sin '0 cos '0 cos^ a sin a cos a sin ■?/' 
sin -0 cos Ip cos^ a cos^ ip cos^ a sin a cos a cos 
sin ip cos a sin a cos ip sin a cos a 



sin^ a 



(52) 
(53) 

(54) 

(55) 

(56) 



F. Astrophysical and cosmological constraints 

Astrophysics has turned out to be quite effective in setting an upper bound on the 077 
vertex in £alp; which therefore holds for the axion as well as for ALPs. 

In the first place, the failure to detect ALPs emitted by the Sun in the CAST experiment 
at CERN has led to 

M > 1.14- 10^°GeV (57) 

for m < 0.02 eV [2Ql. 

On the theoretical side, the most reliable method concerns ALP photo-production through 
the Primakoff process, which takes place when an incoming photon scatters on a charged 
particle and becomes an ALP upon the exchange of a virtual photon. Hot, dense plasmas 
in stellar cores are ideal environments wherein the Primakoff process involving thermal 



16 



photons can occur. Once produced, the ALPs escape because their mean free path is much 
larger than the stellar radius, thereby carrying off energy. Owing to the virial equilibrium, 
the stellar core has a negative specific heat. Therefore it reacts to such an energy loss by 
getting hotter. As a result, the rate of nuclear reactions sharply increases, bringing about 
a substantial change in the observed properties of stars. Since current models of stellar 
evolution are in fairly good agreement with observations, M has to be large enough to 
provide a sufficient suppression of unwanted ALP effects. This argument has been applied 
systematically in a quantitative fashion to the Sun, to main-sequence stars and to red-giants 
stars in globular clusters, with the result [19] 



M > 10^° GeV . (58) 
A consequence of photon-ALP oscillations is that a lower bound on M stronger than 



conditions (57) and (58) - even if much less robust - can be derived for m < 10~^'^eV. In 
this connections, two methods have been put forward. One is based on the observation of a 
time-lag between opposite-polarization modes in pulsar radio emission and yields [53] 

M > 5 ■ 10^° GeV . (59) 

The other involves ALPs emitted by the supernova SN1987A, which would convert them 
into 7-rays in the magnetic field of the Galaxy. Using the absence of these photons in the 
Solar Maximum Mission Gamma-Ray Detector, the lower bound 

ill 



M > 10^^ GeV (60) 



has been derived [54j- We stress however that condition (60 ) is affected by large uncertainties, 
reflecting the lack of precise knowledge of the Galactic magnetic field as well as of the energy 
dependence of the detector response. 

Let us next turn our attention to the cosmological constraints on ALPs. At variance 
with the previous astrophysical analysis, the case of the axion differs drastically from that 
of generic ALPs. 

We recall that cosmology sets strong constraints on the axion properties because of their 
coupling to quarks and gluons (indeed necessary in order to solve the strong CP problem). 
Basically, both thermal and non-thermal mechanisms can produce axions in the early Uni- 
verse. Since this issue is not directly relevant for our discussion, we cursorily summarize the 
main results remarking that the situation is in reality much more complex than sketched 
here. Recalling that fa denotes the scale at which f/(l)pQ is spontaneously broken, only the 
range 

0.6- 10^ GeV < /a < 0.6- 10^^ GeV (61) 

is cosmologically allowed. Thanks to Eqs. ^ and (|6]), this constraint translates into the 
conditions 

10"^ eV < m < 1 eV (62) 

and 

0.7- 10^° GeV < M < 0.7- 10^^ GeV , (63) 



respectively (we have taken for simplicity A; = 1 in Eq. ([6j)). We stress that these bounds 
should be regarded merely as order-of-magnitude estimates. Moreover, the axion is a very 
good candidate for dark matter. More specifically, for m close to 10~^ eV non-thermal 
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production dominates and it behaves as a cold dark matter candidate, whereas for m close 
to 1 eV thermal production dominates and it is a hot dark matter particle [12] . Searches 
for axionic dark matter are currently underway with the ADMX experiment [53]. Finally, 
it has recently been realized that cold dark matter axions ought to form a Bose-Einstein 
condensate [55] . 

Clearly, all these considerations do not apply to ALPs, since they are supposed to interact 
with the rest of the world through the two-photon coupling only. As a consequence, they 
can be produced in the early Universe only thermally through the processes 7 — )■ a 
and e+ e~ — )■ 7 a. It has been shown that in the case m ^ 1 eV - which is the one relevant 
for us as we shall see later - ALPs are relativistic today and their abundance is anyway 
smaller than that of CMB photons [57] . Hence, we are led to the conclusion that the ALP 
considered in this paper are totally unconstrained by cosmology and play no role for the 
dark matter problem. 



III. VERY-HIGH-ENERGY (VHE) BLAZAR SPECTRA 

Among the many achievements of lACTs is the determination of blazar spectra at energies 
above 100 GeV, and to date this task has been accomplished for about 30 sources with 
redshift up to ^ = 0.536 for 3C 279. Most of these blazars are listed in Table |T| 

In view of our later analysis, we carefully address the propagation of a monochromatic 
photon beam emitted by a blazar at redshift z and detected at energy £"0 within the standard 
ACDM cosmological model, so that the emitted energy is Eo{l + z) owing to the cosmic 
expansion. Clearly - regardless of the actual physics responsible for photon propagation - 
the observed and emitted differential photon number fluxes - namely dN/ dE - are related 
by 

$obs(^o, z) = P^^^iEo, z) (^0(1 + z)) , (64) 

where P^^^{Eq, z) is the photon survival probability throughout the whole travel from the 
source to us. We suppose hereafter that Eq lies in the VHE 7-ray band, and throughout 
this Section we employ cgs units for clarity. 



A. Conventional photon propagation 

Within conventional physics the photon survival probability P^]^^{Eq,z) is usually 
parametrized as 

P^%(Eo,^) = e---(^°'^) , (65) 

where t^{Eq, z) is the optical depth, which quantifies the dimming of the source. Note that 
T^{Eo,z) increases with z, since a greater source distance entails a larger probability for 
a photon to disappear from the beam. Apart from atmospheric effects, one typically has 
t^{Eq,z) < 1 for 2; not too large, in which case the Universe is optically thin up to the 
source. But depending on Eq it can happen that t^^Eq, z) > 1, so that at some point the 
Universe becomes optically thick along the line of sight to the source. The value z^ such 



that T^{Eo,Zh) = 1 defines the 7-ray horizon for a given £"0, and it follows from Eq. (65) 



that sources beyond the horizon tend to become progressively invisible as z further increases 



past Zh- Owing to Eq. (65), Eq. (64) becomes 



$obs(^0, Z) = e-^-(^°'^) $em (^o(l + z)) . (66) 
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Source 


Z 


Observed energy range 


Tobs 


Reference Pem 


3C 66B 


0.022 


120 GeV < ^0 < 1.8 TeV 


3.10 + 0.37 


m 


3.00 


Mrk 421 


0.030 


140 GeV < £;o < 6 TeV 


2.33 + 0.22 


m 


2.16 


Mrk 501 


0.034 


150 GeV < £^0 < 6 TeV 


2.09 + 0.20 


m 


1.90 


Mrk 501 


0.034 


150 GeV <Eq< 3 TeV 


2.20 + 0.20 


m 


2.03 


lES 2344+514 


0.044 


180 GeV <Eq< 4 TeV 


2.95 + 0.23 


[62] 


2.70 


Mrk 180 


0.045 


180 GeV < ^0 < 1.4 TeV 


3.30 + 0.70 


[63] 


3.07 


lES 1959+650 


0.047 


190 GeV <Eq< 6 TeV 


2.72 + 0.24 


m 


2.43 


BL Lacertae 


0.069 


170 GeV <Eq< 700 GeV 


3.60 + 0.54 


m 


3.27 


PKS 0548-322 


0.069 


440 GeV < ^0 < 2.2 TeV 


2.80 + 0.32 


[66] 


2.39 


PKS 2005-489 


0.071 


230 GeV < ^0 < 2.3 TeV 


4.00 + 0.41 


m 


3.59 


RGB J0152+017 


0.080 


320 GeV < -Eo < 3 TeV 


2.95 + 0.41 


[68] 


2.47 


W Comae 


0.102 


270 GeV < ^0 < 1.2 TeV 


3.81 + 0.49 


[69] 


3.18 


PKS 2155-304 


0.117 


230 GeV < -Eo < 3 TeV 


3.37 + 0.12 


m 


2.67 


RGB J0710+591 


0.125 


? 


2.80 + 0.30 


m 


? 


H 1426+428 


0.129 


800 GeV < < 10 TeV 


2.60 + 0.61 


m 


0.85 


lES 0806+524 


0.138 


320 GeV <Eq< 630 GeV 


3.60 + 1.04 


m 


2.70 


lES 0229+200 


0.140 


580 GeV < < 12 TeV 


2.50 + 0.21 


m 


0.41 


H 2356-309 


0.165 


220 GeV <Eq< 900 GeV 


3.09 + 0.26 


m 


2.06 


lES 1218+304 


0.182 


180 GeV < ^0 < 1.5 TeV 


3.08 + 0.39 


m 


2.00 


lES 1101-232 


0.186 


280 GeV < ^0 < 3.2 TeV 


2.94 + 0.20 


m 


1.72 


lES 0347-121 


0.188 


300 GeV < ^0 < 3.0 TeV 


3.10 + 0.25 


m 


1.87 


lES 1011+496 


0.212 


160 GeV < ^0 < 600 GeV 


4.00 + 0.54 


m 


2.90 


S5 0716+714 


0.31 


180 GeV <Eq< 680 GeV 


3.45 + 0.58 


m 


1.60 


PG 1553+113 


0.40 


95 GeV < £^0 < 620 GeV 


4.27 + 0.14 


m 


2.48 


PKS 1222+21 


0.432 


80 GeV < £o < 360 GeV 


3.75 + 0.34 


[82] 


2.47 


3C 66A 


0.444 


230 GeV <Eq< 470 GeV 


4.10 + 0.72 


[83] 


1.28 


PKS 1424+240 


0.5 


140 GeV < ^0 < 500 GeV 


3.80 + 0.58 


m 


1.16 


3C 279 


0.536 


80 GeV < £o < 480 GeV 


4.10 + 0.73 


[85] 


2.05 


APem 










3.18 


(Pern) 










2.22 



TABLE I: Blazars observed so far with the lACTs with known redshift z, measured energy range, 
measured spectral index Pobs) ^-nd unfolded spectral index at emission Pem using the FRV model 
of the EBL. Statistical and systematic errors are added in quadrature to produce the total error 
reported on the measured spectral index. When only statistical errors are quoted, systematic errors 
are taken to be 0.1 for H.E.S.S. and 0.2 for MAGIC. The last two rows show the spread APgm of 
the values of Pem discarding errors and the average value (Pem) including errors, respectively. 



Whenever dust effects can be neglected, photon depletion arises solely when hard beam 
photons of energy E scatter off soft background photons of energy e permeating the Universe 
and produce e^e~ pairs through the standard 77 — > t^t~ process. Needless to say, in order 
for this process to take place enough energy has to be available in the centre-of-mass frame 
to create an e^e~ pair. Regarding E as an independent variable, the process is kinematically 
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allowed for 



2 ml 



(67) 



E {1 — cos if) ' 

where denotes the scattering angle and me is the electron mass. Note that E and e change 
along the beam in proportion of I + 2;. The corresponding Breit- Wheeler cross-section is 



e, ^) - 1.25 ■ lO'^^ (l - (3^) 



2/3 2) + (3-/3^) ln(i±|) 



cm 



(68) 



which depends on E, e and (p only through the dimensionless parameter 



/3{E,e,^) = 



1 - 



2 mid 



2 ^4 



Ee (1 — cosy^) 



1/2 



(69) 



and the process is kinematically allowed for > 0. The cross-section a.y.y{E, e, ip) reaches its 
maximum a!^^^ ^ 1.70 ■ 10~^^ cm^ for P ~ 0.70. Assuming head-on collisions for definiteness 
(93 = tt), it follows that a^^{E, e, tt) gets maximized for the background photon energy 



e{E) 



500 GeV 
E 



eV 



(70) 



where E and e correspond to the same redshift. 

Within the standard ACDM cosmological model t^{Eq,z) arises by first convolving 
the spectral number density n^{e{z),z) of background photons at a generic redshift with 
a^^{E{z), e{z), ip) along the line of sight for fixed values of z, ip and e{z), and next integrat- 
ing over all these variables [32]. Hence, we have 



Ty{Eo, z) 



dz 



dl{z) 



d(cos p) 



1 — cos (p 



X 



(71) 



de{z) n^{e{z), z) a^^{E{z),e{z), ip) , 

where the distance travelled by a photon per unit redshift at redshift z is given by 

dl{z) c 1 



dz 



(1 + z) [VtK + VLm (1 + zY] 



31 1/2 



(72) 



with Hubble constant Hq ~ 70 Kms~^ Mpc~^, while — 0.7 and VLm — 0.3 represent 
the average cosmic density of matter and dark energy, respectively, in units of the critical 
density pcr — 0.97 ■ lO^^^gcm"^. 

Once n^{e{z), z) is known, t^{Eq, z) can be computed exactly, even though in general the 
integration over e{z) in Eq. (71) can only be performed numerically. 

Finally, in order to get an intuitive insight into the physical situation under consideration 
it may be useful to discard cosmological effects (which evidently makes sense for z small 
enough). Accordingly, z is best expressed in terms of the source distance D = cz/Hq and 
the optical depth becomes 

^^(^•^) = a;M' '''' 
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where Xj{E) is the photon mean free path for 77 — )■ e~^e referring to the present cosmic 
epoch. As a consequence, Eq. (65) becomes 



and so Eq. (66) reduces to 



Note that we have dropped the subscript for simphcity. 



(74) 



(75) 



B. Extragalactic Background Light (EBL) 

Blazars detected or detectable in the near future with lACTs he in the VHE range 



100 GeV < Eq < 100 TeV, and so from Eq. (70) it follows that the resulting dimming is 



expected to be maximal for a background photon energy in the range 0.005 eV < eo < 5eV 
(corresponding to the frequency range 1.21 ■ 10^ GHz < z/q < 1.21 ■ 10^ GHz and to the 
wavelength range 2.48 /im < Aq < 2.48 ■ 10^ /im), extending from the ultraviolet to the far- 
infrared. This is just the EBL. We stress that at variance with the case of the CMB, the 
EBL has nothing to do with the Big Bang. It is instead produced by stars in galaxies during 
the whole history of the Universe and possibly by a first generation of stars formed before 
galaxies were assembled. Therefore, a lower limit to the EBL level can be derived from 
integrated galaxy counts [86] . 

Determining the spectral number density n^{e{z),z) of the EBL is a very difficult task. 
It is affected by large uncertainties, arising mainly from foreground contamination produced 
by zodiacal light which is various orders of magnitude larger than the EBL itself [87J. While 
it is beyond the scope of this paper to review the various attempts in this direction, we 
briefly summarize below the present situation. Basically, six different approaches have been 
pursued: 

• Forward evolution - This is the most ambitious approach, since it begins from first 
principles, namely from semi-analytic models of galaxy formation in order to predict 
the time evolution of the galaxy luminosity function 



Backward evolution - This starts from observations of the present galaxy population 
and extrapolates the galaxy luminosity function backward in time. Among others, 
this strategy has been followed by Stecker, Malkan and Scully (SMS) [90j and by 
Franceschini, Rodighiero and Vaccari (FRV) 



Inferred evolution - This models the EBL by using quantities like the star formation 
rate, the initial mass function and the dust extinction as inferred from observations [39l 

m 



Minimal EBL model - This relies upon the same strategy underlying the previous item 
but with the parameters tuned in order to reproduce the EBL lower limits from galaxy 
counts [92]. 

Observed evolution - This method has the advantage to rely only upon observations by 
using a very rich sample of galaxies extending over the redshift range < z < 1 [93j . 
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• Compared observations - This technique has been implemented in two ways. One 
consists in comparing observations of the EBL itself with blazar observations with 
lACTs and deducing the EBL level from the VHE photon dimming [94J . Another starts 
from some 7-ray observations of a given blazar below 100 GeV where EBL absorption 
is negligible and infers the EBL level by comparing the lACT observations of the same 
blazar with the source spectrum as extrapolated from the former observations [95] . In 
the latter case the main assumption is that the emission mechanism is presumed to 
be known with great accuracy. In either case, the crucial unstated assumption is that 
photon propagation in the VHE band is governed by conventional physics. 

As it is evident, the latter approach does not apply to the DARMA scenario, and so it will 
not be considered any further. As far as the backward evolution approach is concerned, the 
models of SMS predict a much higher EBL level as compared to the model of FRV. Recently, 
the SMS models have been ruled out by Fermi/LAT observations [96]. On the other hand, a 
remarkable agreement exists among the FRV model and the other models based on forward 
evolution, inferred evolution and observed evolution. Throughout this paper, we adopt the 
FRV model mainly because it supplies a very detailed numerical evaluation of the optical 



depth based on Eq. (71), which will henceforth be denoted by t^^^{Eq, z) [97]. Regretfully, 
the errors affecting rJ^^^(i?o, z) are unknown. 



C. Understanding observed VHE blazar spectra 

As a preliminary step to find out the potential relevance of the DARMA scenario for 
available observations of VHE blazars, we consider in some detail the energy range 0.2 TeV < 
Eq < 2 TeV where most of the blazars in question have been detected. It follows from Eq. 



(70) that the EBL energy band where cr^^(-E'o, eo, vr) becomes maximal is 0.25 eV < eo < 
2.5 eV (corresponding to 6.07- 10^ GHz < uq < 6.07-10^ GHz and 0.50 /im < Aq < 4.94 /xm). 

So far, two specific processes have been proposed which can give rise to the emission of 
VHE photons [S3]. 

• One is the synchro- s elf- Compton (SSC) mechanism, according to which relativistic 
electrons first emit X-ray photons by spiralling in the source magnetic field, which are 
subsequently boosted into the VHE 7-ray band by inverse Compton scattering off the 
parent electrons. In some cases, also external electrons can substantially contribute to 
the inverse Compton. 

• The competing mechanism is hadronic pion production (HPP) in proton-proton scat- 
tering, with neutral pions immediately decaying into VHE 7-ray pairs. 

It turns out that both mechanisms lead to emission spectra which are so far observation- 
ally indistinguishable, and in particular within the energy range 0.2 TeV < E < 2 TeV they 
both predict a single power-law behaviour for blazar emitted spectra 

^,^{E) = K E-^^- , (76) 

where i^' is a suitable constant. 

We next turn our attention to the observed energy spectra. It follows directly from 
observations that blazar spectra are successfully fitted by a single power law [98j 

<l>obs(^o,^)=^^o"^°"^'^ ■ (77) 
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FIG. 2: The observed values of the observed spectral index Fobs versus the source redshift for all 
blazars detected so far in the VHE band are represented by dots and corresponding error bars. 



As a consequence, the observed spectra of all blazars detected so far are characterized by 
the observed spectral index Fobs, which is reported in Table [T] for every source. It is also 
very useful to plot Fobs versus the source redshift z for all detected VHE blazars in Figure 
[2] (blobs with error bars). 

Let us first try to understand what Figure |2] is telling us leaving aside any theoretical 
prejudice. A striking feature is that the horizontal strip 3.5 < Fobs < 4.5 is almost uniformly 
populated for all considered redshifts, which would suggest that Fobs is independent of z. 
However, things are different for the lower strip 2.5 < Fobs < 3.5. Because it is populated 
only up to z ^ 0.2 (with the exception of a single source close to z ~ 0.3), the above inter- 
pretations is ruled out and we are forced to conclude that Fobs correlates with z. Actually, 
when looking at Figure [2] from this viewpoint a simple trend is easily recognized: Fobs in- 
creases linearly from 2.5 - 3 at 0.1 < 2; < 0.2 to roughly 3.5 - 4 at 0.3 < z < 0.6. Similarly, 
also for Fobs > 3 a linear increase is found - even if with a different slope - but the number 
of sources with Fobs > 3 following this behaviour decreases as z increases until it vanishes 
for z > 0.25. 

A qualitative understanding of this situation emerges naturally by taking the EBL at- 
tenuation into account. We stress in the first place that rather nearby blazars - such as 
those at z < 0.05 - do not practically suffer EBL absorption at the energies probed so far, 
thereby implying that the shape of their observed VHE spectra should be the same as that 
of the emitted spectra, namely Fobs — Lem- This is an important fact, since it allows us to 
see directly the blazar spectra at emission. In addition, we show in Appendix B that an 
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approximate analytic expression for the optical depth within the FRV model is given by 

T^^^iEo, z) ^ 2.25 a ( ^^^^ ) I{z) , (78) 



500 GeV 



with 0.9 < a < 3.6 and roughly I[z) ~ z. Hence, by combining Eqs. (66), (76) and (78) the 
expected observed flux is 



«(^o,^) = i^exp|-2.25a(^^^^y'''/(z)| E,^-- (1 + z)" 



Now, Eq. (79) possesses to two conceptually distinct implications: 



^obs (-^0; z) is exponentially damped as the energy increases, thereby entailing that it 
gets much softer than the emitted flux. 

^^^{Eq, z) is exponentially suppressed as the distance increases, so that sufficiently 
far-away sources tend to become invisible. 



Although Eq. (|79|) holds up to 2 TeV only, one can check that these conclusions remain true 



under the replacement r^PP(£'o,2;) — )■ {Eq, z) up to 100 TeV. 



We can relate the expected observed spectral index E^^^^^) to Egm and z by best-fitting 



the l.h.s. of Eq. (79) to the power-law expression (77) over the energy range where the 
considered source is observed, which is reported in Table [1} 

Because I{z) is independent of Eq, it is unaffected by the considered best-fitting proce- 
dure, and so we have r|^^g(2;) ~ I[z) ~ z up to logarithmic corrections. This indeed explains 
in a qualitative fashion both why nearby sources with 2.5 < Fobs < 3 get replaced by sources 
with 3.5 < Fobs < 4 at larger redshift according to a linear trend and why nearby sources 
with Fobs > 3 follow a similar linear behaviour up to a point where Fobs would be so large 
that the source becomes invisible at sufficiently large distances, thereby disappearing from 
Figure |2] for z large enough. 

An intrinsic correlation - due to an observational bias - between the spectral index and 
the distance, however, cannot be excluded on the basis of the experimental data. In addition, 
given the blazar sequence [99j, the fact that the Inverse Compton bump moves according 
to luminosity might give a bias related to the fact that for the same energy range we are 
actually sampling different regions of the spectral energy distribution. Finally, the fact that 
the upper limit of the energy sampled decreases with energy might introduce in itself a bias. 
However, a direct search for spectral index hardening associated with blazar variability gave 
no evidence |100j . Attempts are presently done (see for example [lOlj ) to analyze individually 
blazars and derive their spectral energy distribution from multi-wavelength data. We hope 
that in a near future this work will be made more precise. However, the models have 
presently a large uncertainty, since, in order to have reasonably constrained fits, one must 
assume a purely leptonic emission and a 1-zone SSC emission mechanism, while we have 
indications that the situation can be more complicated for most blazars we know in detail. 

In order to derive the exact value of Fem for the various blazars from observations the 



use of Eq. (79) with <^^{Eq, z) — (^ohs{Eo,z) combined with Eq. (77) would be unsuited 
because of its approximate character. A better strategy consists in first de-absorbing Fobs 



for every source by employing Eq. (66) with t^{Eq,z) — )■ t^^[Eq,z) combined with Eq, 



(77), and next inferring Fem by best-fittmg the resulting $em(-S) to the power-law expression 
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(76) over the energy range where the considered source is observed (see Table [T]). Observe 
that since Fem depends hnearly on Fobs, the derived values of Fem have the same error bars 
of Fobs as reported in Table [l] to the extent that errors in t^^"^ {Eo, z) are neglected (they 
are actually unknown because they are not quoted by the authors). We are of course well 
aware that the correct procedure would be to first de-absorb each point of the observed 
spectrum of a given source and next best-fit these points to a power-law. Unfortunately, 
the observed energy points with related error bars are not available from published papers, 
and this explains why we have simply de-absorbed Fobs - hence Fem has to be understood as 
the average emitted spectral index for the source in question - but we want to remark that 
for our purposes this simplified approach is adequate. The same strategy has been used for 
a different model of the EBL [89]. Our results are listed in Table [l| along with the spread 
AFcm of the value of Fgm discarding errors and the average value (Fem) including errors. 



D. The cosmic opacity problem 

It is evident from Table [Tjthat the values of Fem for some far-away sources are considerably 
smaller than those for nearby blazars, with the exception of the two sources H 1426-1-428 
and lES 0229+200. 

As a consequence, the cosmic opacity problem arises concerning the physical mechanism 
responsible for such a behaviour involving Fem and z. 

Certainly cosmology does not help, because no important evolutionary effect is expected 
to take place for redshifts up to ^ ~ 0.54 at which the most distant blazar has been detected. 

Alternatively, one might guess that it is due to a volume selection effect, since intrinsically 
brighter sources are the exception rather than the rule. However, the emitted flux depends 
not only on the slope but also on the normalization factor, which varies by three orders of 
magnitude over the sample of considered sources. Actually, the existence of the two rather 
nearby blazars H 1426+428 and lES 0229+200 with the hardest emitted spectrum explicitly 
shows that small Fem does not mean large z. So, also this attempt is unsatisfactory. 

Yet another possible explanation consists in assuming that far-away blazars are intrinsi- 
cally different from nearby ones, but to the best of our knowledge no convincing explanation 
of this circumstance has been put forward so far. 

One might also argue that a solution could come from the fact that some observed blazars 
are in a quiescent state whereas others are flaring. In fact, because of EBL absorption we 
might be seeing progressively more distant blazars only during stronger flares (but not all 
distant blazars are flaring, like e.g. PG 1553+113 which has been observed to have nearly 
the same luminosity for five years). As a consequence - working within the SSC mechanism 
for definiteness - we could run the risk to compare Fem for different sources at different 
positions on the Compton peak, since flaring causes this peak which normally lies below 
100 GeV to slightly shift towards higher energies |102] . Clearly, the slope near the bottom 
of the pick is steeper than close to the tip, and this circumstance would produce a harder 
emission spectrum for flaring sources. However, such a possibility seems to us quite unlikely. 
For, the observed energy range of flaring sources is generally considerably wider that the 
width of the peak |102] and observations above 100 GeV invariably show that a single power 
law behaviour provides an excellent fit to the data. Hence, we see that we are inferring the 
spectral index well below the pick whether or not a flare takes place. 

In conclusion, no satisfactory explanation for the considered behaviour involving Fem and 
z seems to emerge. 
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As a matter of fact, shifting from the astronomical to the physical point of view makes the 
issue more clear-cut. It is known since a long time that Fem = 1-5 arises from the first-order 
Fermi acceleration mechanism with newtonian shocks for an electron injection spectrum 
equal to 2 |1U3] . For this reason, when the cosmic opacity problem was first perceived in 
2006 it was thought that the inferred values of Fgm were too small to agree with conventional 
physics assuming current EBL models [36], and indeed values Fem < 1-5 were considered 
unphysical e.g. by the H.E.S.S. collaboration. However, it has recently been shown that the 
required low values of Fem can be achieved in the presence of strong relativistic shocks jlQ] , 
because of photon self-absorption inside the source jTlJ or by the inverse Compton scattering 
of CMB photons by shock-accelerated electrons in the jet |12]. 

While these results are gratifying, one still wonders why this kind of physical effects are 
important for distant blazars only. Just as before, an answer seems hardly in sight. 

So, either way it is argued no simple solution to the cosmic opacity problem emerges 
within conventional physics. 

Now, what the discussion in Subsection III-C has shown is that the origin of this cosmic 
opacity problem is not qualitative - the 2;-dependence of Fobs(-2) comes out right - but purely 
numerical, namely because the EBL level predicted by conventional physics is too high. Were 
the EBL level somewhat smaller, the cosmic opacity problem would automatically disappear. 

A way out of the cosmic opacity problem appears therefore to call for some sort of 
unconventional photon propagation which ultimately reduces the cosmic opacity arising 
from the EBL, thereby bringing the values of Fgm for distant blazars in the same ballpark 
of those for close ones. 

Various attempts at reducing the cosmic opacity along these lines have been proposed 
and they are schematically summarized below: 

• A revolutionary option contemplates a breakdown of Lorentz invariance jl04] . 

• An alternative possibility concerns the emission of cosmic rays from blazars - rather 
than photons - with energy smaller than 50 EeV. These cosmic rays can travel unim- 
peded over cosmological distances and they can interact with the EBL well before 
reaching our galaxy. In such an interaction secondary photons are produced, that are 
ultimately detected by the lACTs |105j . 

• A different proposal relies upon photon-ALP oscillations, which requires the presence 
of magnetic fields somewhere along the line of sight. As already pointed out, two 
concrete realizations of this idea have been investigated. One of them - the DARMA 
scenario p^H26] - assumes that photon-ALP oscillations take place during propagation 
in intergalactic space, where large-scale magnetic fields in the nano-Gauss range are 
supposed to exist. Large-scale magnetic fields of this strength are consistent with 
current upper bounds and even with the results of the AUGER observatory (more 
about this, later). The other is in a sense complementary, because it presupposes a 
7 — )■ a conversion inside the blazar and a a — )■ 7 conversion in the Milky way [27] . 
Although the properties of the Galactic magnetic field are rather well known, those of 
the magnetic field in the blazar are not, and so it is not clear whether the first step of 
this mechanism actually takes place and if so how large is its efficiency [29] . 




The present paper is devoted to a careful investigation of the solution based on the 
DARMA scenario, even if its scope is by far more general. 
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IV. DARMA SCENARIO 



Our aim is to offer a detailed description of tlie structure of tlie DARMA scenario, and 
in particular to show how the photon survival probability P^_^^^^{Eo, z) can be computed 
in terms of the properties of the intergalactic medium in which the photon/ ALP beam 
propagates. 



A. An intuitive insight 

We find it instructive to restate in a slightly different fashion the reason why the mech- 
anism of photon-ALP oscillations allows to substantially reduce the EBL absorption. We 
neglect here cosmological effects for simplicity. 

We suppose that VHE photons are both emitted and detected as usual, but that along 
their way to us they convert into ALPs and back into photons. Accordingly, the number 
Nc of either 7 — )■ a or a — )■ 7 conversions must necessarily be even, and we may schemat- 
ically regard the beam propagation in large-scale magnetic fields as a succession of such 
conversions. Assuming ideally that each conversion occurs suddenly at some space point, 
the source distance D gets divided into a number iVc + 1 of steps of equal length L, over 
which a beam particle behaves either as a real photon or as an ALP. Hence, a beam particle 
exhibits an overall behaviour as a real photon over a total length equal to 

^ 2(JV, + 1) V 2N, + 2j ^ ' 

We intuitively expect Nc to increase with the photon-ALP oscillation probability - and so 
with B/M - which leads in turn to a slight decrease of starting from D. Correspondingly, 



since ALPs do not suffer EBL absorption Eq. (74) gets presently replaced by 



P^_^^^^{E,D) = e-^^/^^(^) = e"(^"^^)^/^^(^) , (81) 

and thanks to its exponential dependence on D^/X^{E) even a small decrease of start- 
ing from D produces a large enhancement of PP^^^^{E,D) as compared to PS^ {E,D) 



referring to conventional physics and given by Eq. (74) |lU6j . 



B. General strategy 

Our ultimate goal consists in the evaluation of the photon survival probability 
Plf^^^^^{E(), z) from a blazar at redshift z to us when allowance is made for photon-ALP 
oscillations as well as for EBL photon absorption. It is indeed clear that the considerations 
developed in Section III can be extended to account for photon-ALP oscillations by the 
replacement P^X^{Eo,z) ^ P^^^^'^iEo, z). 

An exact treatment would however be impossible because of the large uncertainty affect- 
ing the configuration of the magnetic field ultimately responsible for photon-ALP oscilla- 
tions. 

As a matter of fact, the line of sight to a distant blazar is expected to traverse magnetic 
fields extending over a variety of scales. A magnetic field is certainly present inside the 
source |107] . Furthermore, the Milky Way magnetic field can give a nontrivial contribution 
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to the effect under consideration |108j . and the same is true if the hne of sight happens to 
cross a cluster of galaxies because intracluster magnetic fields are known to exist with a 
strength similar to that of the Galactic field jl09| . Finally, large-scale magnetic fields can 
play a key role fTTU| llllj . Here, our attention will be restricted to magnetic fields of the 
latter sort. 

Unfortunately, almost nothing is known about the morphology of large-scale magnetic 
fields, which reflects both their origin and the evolutionary history of baryonic matter. While 
it is evident that their coherence length cannot be arbitrarily large, no reliable estimate of 
its value is presently available. As far as our analysis is concerned, this means that we 
cannot suppose that large-scale magnetic fields are homogeneous over the whole distance 
to the source, but their spatial dependence is largely unknown. The usual way out of this 
difficulty amounts to suppose that large-scale magnetic fields B have a domain-like structure. 
That is, B is assumed to be homogeneous over a domain of size Ldom equal to its coherence 
length, with B randomly changing its direction from one domain to another but keeping 
approximately the same strength |110l 1111] . 

Thus, the whole propagation process of the considered photon/ALP beam can be recov- 
ered by iterating the propagation over a single domain as many times as the number of 
domains crossed by the beam, taking each time a random value for the angle ip between B 
and a fixed fiducial direction equal for all domains. In this way, we are effectively led to the 
much easier problem of photon-ALP oscillations in a homogeneous magnetic field (solved 
exactly in Subsection II-E). 

What still remains to be done at this point is to take photon absorption into account. 
This is easy because photon absorption is independent of the properties of the photon-ALP 
oscillation mechanism and vice-versa. 

So, our strategy can be implemented through the following steps [21j: 

• We work within the strong-mixing regime so as to ensure that the photon-ALP oscil- 
lation probability is both maximal and energy- independent. We have seen that such 



a condition requires E > E^, with the energy threshold E^ defined by Eq. (44). But 
demanding the strong-mixing regime to take place for E > 100 GeV evidently requires 
E^ < 100 GeV, which sets an upper bound on the ALP mass. 

• We evaluate the transfer matrix across the generic n-th domain W„(-E'o) V'n); where 
ipn accounts for the random orientation of B in the domain in question. Note that 
Un{EQ,ipn) depends on Eq only because of the energy-dependence of EBL absorption. 

• Iteration of the latter result over the total number of domains crossed by the 
beam from the blazar to us yields the total transfer matrix V({Eq, z; ipi, from 
which the photon survival probability P^^^{Eq, z;ipi, ...,ipNa) can be computed for 
fixed values of the angles i/ji, ...,1^^^ every domain. 

• Finally, P}^_^^^^{Eq, z) emerges by averaging P.y^^{EQ,z;ipi, ...jipiyj over all angles 

Our discussion is framed within the ACDM cosmological setting, and so the redshift z 
is the obvious parameter to express distances. Because the proper length per unit redshift 



at redshift z is still given by Eq. (72), a generic proper length extending over the redshift 
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interval [za, Zb] {za < Zb) is 

, dl(z) , n dz , , 

L{za, Zb) = / dz ~ 4.29 ■ 10^ / — Mpc , (82) 



' A. (1 + ^) [0.7 + 0.3(1 + 2)']'^' 

which approximately reads 

(1 + 1 45 z \ 
, , ' Mpc . (83) 

1 + 1.45 2;^/ 

This result will be applied in particular to evaluate the size of the magnetic domains. 

C. Photon absorption 

We proceed to extend the discussion in Subsection II-E so as to take EBL absorption into 
account. This task is greatly facilitated by the fact that the latter effect is independent of 
the photon-ALP conversion mechanism. 

We have seen that the propagation of a monochromatic photon/ ALP beam is formally 
described as a three-level non-relativistic quantum system with Hamiltonian Hq given by 



Eq. (21 ) and expressed in terms of the mixing matrix JHq. Taking advantage from this fact, 
the inclusion of EBL absorption amounts to suppose that the photon/ ALP beam is actually 
analogous to an unstable quantum system with decay probability 

Pdecay = e-^/^^^^^ , (84) 

where X'y{E) denotes the photon mean free path. As is well known, such a decay probability 
arises from the inclusion of an absorbitive term — Aabs into the Hamiltonian, with 

^^^^ ^ 2XjE) ■ ^^^^ 



More specifically, since photons undergo absorption but ALPs do not, ^ in Eq. (27) 



becomes 

/Api + Aabs \ 

-M^"^ = Api + Aabs A,^ , (86) 

V A,^ AaaJ 

where - in parallel with the treatment of Subsection II-E - we are first supposing that B 
lies along the z-axis. 

As we said, we work throughout within the strong-mixing regime and therefore condition 
43 ) has to be met. Recalling the explicit expression for the various A-terms entering Eq. 
86) and defined in Subsection II-E, A^^^^ takes the simpler form 



Aabs 










Aabs 







A(j^ 


/ 



= Aabs A,^ , (87) 



which is denoted by the same symbol for notational simplicity (only Eq. (87) will be used 
hereafter). Note that m and Upi presently drop out of A^*^°\ Just as before, use of the results 
contained in Appendix A with Ai ^ A4^^^ directly gives the corresponding eigenvalues 
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A2 
A3 



4A^(E) 



4A^(E) 



(^1-^1-452) 
(^1 + ^1-452) 



where we have set 



5 = 



BX^{E) 
M 



(89) 
(90) 

(91) 



roughly measuring the ratio of the photon mean free path to the photon-ALP oscillation 
length. Hence, the transfer matrix associated with the reduced Schodinger-like equation 
m\ with Mq reads 



U{y, yo; 0) = e'^'^^^-^"'^ Ti(0) + e'^'^y-^'''^ T2(0) + e^^^^^"^") r3(0) 



(92) 



with the matrices Ti(0), T2(0) and T^iO) dictated by Eqs. (166), (167) and (168) as special- 
ized to the present case. Explicitly 
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(93) 



(94) 



(95) 



Clearly, we will need the generalization of this result to the case in which B forms an 
arbitrary angle ip with the ^-axis. Proceeding exactly like in Subsection II-E, we find 



where the matrix V{ip) is given by Eq. (50). This yields 

M = 



^abs 







^abs 



Aay sinip 

COSIp 





sinip Aay cosip 
and now the resulting transfer matrix evidently reads 

U{y,yo;ij) = V\i;)Uiy,yo;0)V{^P) 



whose explicit form arises by inserting Eq. (92) into Eq. (98). We obtain 
U{y, yo; ^) = e^^^^^"^") Ti{^) + e'^^^^-y^^ T2{tp) + e'^^^y-^''^ 

with 

cos^ tp — sin ip cos -ip 
Ti (■?/') = I — sin ip cos %p sin^ %p 




(96) 

(97) 

(98) 
(99) 
(100) 



30 



2 Vl-4<5 2 
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sin cos 
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(101) 



(102) 



Vi-452 

We stress that due to the imaginary nature of A^bs the mixing matrix Ai is not self- 
adjoint, and so the transfer matrix U{y,yo]ip) fails to be unitary. In addition, the Von 



Neumann-like equation (46) becomes 



= pM^ - Mp . 

dy 



(103) 



Still, it is straigthforward to check that Eq. (47) remains valid with Uo{y,yo) — >• W(|/,?/o; V') 
in spite of the fact that A^^ 7^ Ai. Since we now have p^yY 7^ p{y): it follows that in general 
TTp{y) 7^ 1, and so the probability that a photon/ALP beam initially in the state pi will be 
found in the state p2 after a distance y is presently given by 



P,,^,, (y) = Tr (p2 Uiy, 0; V^) pi (y, 0; ?A)) 



(104) 



where we assume Trpi = Trp2 = 1 as before. 



D. Intergalactic medium (IGM) 

As is well known, the absence of the Gunn-Peterson effect [112] is generally regarded 
as evidence that the IGM is ionized, and from the resulting high electrical conductivity it 
follows that the electron number density ne{z) traces the cosmic mass distribution. Because 
of this fact, we have 

ne{z) =ne,o{l + S{z)){l + zf , (105) 
where S(z) = {p(z) — p{z)) / p{z) is the mass density contrast and fie.o is the average electron 



number density. As a consequence, Eq. (25) entails for the plasma frequency 



^pi(^) = c^pi,o(l + S{z)f' (1 + zf , (106) 

with ojpifi obviously corresponding to fiefl. 

Observations of the primordial abundance of the light elements yields n^fi — 1.8 • 
lO^^cm"^, but it has been argued that in the z < 1 Universe which is relevant for us 



ne{z) ought to be smaller than Ugfl by a factor 15 |113] . Correspondingly, from Eqs. (25) 



and (106) we get 

Upi{z) ~ 4.04 ■ 10^^^ (1 + zf^^ eV , (107) 

where the (1 + 5{z)Y^'^ factor has been dropped because irrelevant. 

A crucial issue concerns the large-scale magnetic fields traversed by the beam, whose 
origin and structure is still unknown to a large extent. A possibility is that very small 
magnetic fields present in the early Universe were subsequently amplified by the process of 
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structure formation |114] . An alternative option is that the considered magnetic fields have 
been generated in the low-redshift Universe by energetic quasar outfiows |115j . Finally, it 
has been suggested that large-scale magnetic fields originated from the so-called Biermann 
battery effect |116] , namely from electric currents driven by merger shocks during the struc- 
ture formation processes. Presumably, all these mechanisms can take place, even if it is 
presently impossible to assess their relative importance |110l 1111] . At any rate, we suppose 
that magnetic fields already exist out to the redshift z = 1. 

Owing to the high conductivity of the IGM, the magnetic fiux lines can be thought as 
frozen inside the IGM. Therefore, fiux conservation during the cosmic expansion entails that 
B scales like the volume to the power 2/3, thereby implying the magnetic field strength in 
a domain at redshift z is |lllj 

B = Bo {l + zf . (108) 

In 2007 the AUGER collaboration reported positive evidence for a correlation between 
charged cosmic rays and candidate sites for emission |117j . More in detail, the AUGER 
collaboration found that 20 out of the 27 recorded events with energy larger than 57EeV 
are located within 3.1° of an AGN closer than 75Mpc from Earth. The conclusion drown 
from the AUGER collaboration is that such a result is inconsistent with the hypothesis of 
an isotropic distribution of these cosmic rays with at least a 99 % confidence level from a 
prescribed a-priori test. As explained elsewhere, this fact supports the existence of large-scale 
magnetic fields with coherence length Ldom in the range 1 — 10 Mpc and strength Bo in the 
range 0.1 — 1 nG ai z = |118j . However, such a correlation has become considerably weaker 
when a larger data set (69 events, including the events on which the previous publication was 
based) has been recorded and analyzed by the AUGER colaboration |119j . The fraction of 
events correlated to a nearby AGN is in the most recent publication of 38 % - to be compared 
with an expected value of 21 % in the case of no correlation - and no a-priori probability 
estimate is provided in the new paper. The present situation appears to us unclear, even 
though the conclusions obtained from the first AUGER results are still statistically consistent 
with the more recent ones. 

For this reason, we prefer to avoid committing ourselves with any conclusion relying 
upon the AUGER data and we consider only well-established upper bounds. They depend 
on the size of their domain-like structure Ldom and within the current cosmological setting 
they take the form |120j 

5o < 3.8nG for Ldom = 50Mpc , (109) 



fio < 6.3nG for Ldom = IMpc . (110) 
It is usually supposed that 1 Mpc < Ldom < 10 Mpc, and so we will assume throughout 

5o<6nG. (Ill) 

Within the cosmological context, the overall structure of the cellular configuration of 
large-scale magnetic fields is naturally described by a uniform mesh in redshift space with 
elementary step Az, which can be constructed as follows. The magnetic domain closest to 
us and labelled by n = 1 extends from to Az. Hence, its size i^dom '^^^ written as 

(1) 



41 = ^(0, A^) = I ^ I 5 Mpc , (112) 
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where L(0, Az) is the domain's proper size a nd th e second equahty stresses our preferred 
choice for L^^^^- By combining Eqs. (83) and (112) we infer 



Az ~ 1.17-10" 



r{l) 

5Mpc 



(113) 



indeed in agreement with the hnear Hubble law. Because our mesh in redshift space is 
uniform, Az sets the redshift size of all magnetic domains. Therefore for a source at redshift 
z the total number of magnetic domains crossed by the beam can be estimated as 



z 



d- — - 0.85 ■ 10^ 



Az 



5Mpc 

^dom 



(114) 



and since we are assuming z < 0.54 we have < 0.46 ■ 10*^(5 Mpc/L[^^^). Furthermore, the 
n-th domain extends from z = [n — l)Az to z = nAz and its proper size can be written as 



^dom 



L((n — l)Az,nAz). Thanks again to Eq. (83), it reads 



-^dom — 



2.96 • lOMn 1 + 



1.45 Az 



1 + 1.45 (n - l)Az 
4.29 ■ lO^Az 



Mpc 



(115) 



Mpc 



1 + 1.45 {n - l)Az 

where the last equality is justified by the fact that our analysis is confined to 2; < 0.54 in 



conjunction with Eq. (114). 



E. Propagation over a single domain 

We are now in position to describe the propagation of the considered photon/ALP beam 
across the n-th magnetic domain. 

The transfer matrix is directly supplied by Eq. (99), which for notational convenience we 
rewrite as 



(116) 



with 



T2(^„) = 



( -l+v/l-4,5^ 
-1+^1-4^2 



COS^ — sin ll)n COS 

— sin COS 'j/'n siu^ "^ri 





sin^ t/;^ 



sin cos '?/'„ 



sin-?/' 

n 

COS 



^ COS^ 7/'„ 



2^1^ 



sin^/'n 



a/i-45„: 



COS '^ri 

















sin " 


v/l-4'5^ 


iSn 


cos 

-4<52 

n 




2^1- 


-452 



;ii7) 



(118) 
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■ -z I l+\/l-4<52 . . \ 

/ Sm , Sm iPn COS , " Sill Ipn 



sm ihn COS ibn — \ COS^ ibn —Ms^ 



COS Ipr, 



sin 



cos 



-1+^1-4^2 



(119) 



where ipn denotes the angle between and the z-axis, which is fixed for all domains. 
Moreover, we have set 



X 



in) 



A 



(n) 



4A^ 



(l-^r^45l) 



with 



En = Eq 



5n = 



1 + (n - 1) A^) 

Bn\^;\Eo) 
M 



and we have introduced the shorthand 

A(")(Eo) = A^ (En) 



In addition, Eq. (108) implies 



Bn = Bn 



l + {n-l)Az) 



(120) 
(121) 
(122) 

(123) 
(124) 

(125) 
(126) 



What remains to be done in to evaluate the photon mean free path X^'y^\EQ). A convenient 
procedure is as follows. Let us suppose to observe two hypothetical sources located at both 
edges of the ra-th domain. Then we apply Eq. (66) to either source. With the notational 
simplifications $obs(-Eo, z) $(-E'o) and $em {Eo{l + z)) $ {Eo{l + z)), we have 

$(Eo) = e-"^(^«'("-i)^^) $ (En) , (127) 

$(^0) = e-"-(^0'"^") $ (E„,+i) , (128) 
and so the fiux change across the considered domain is 

$ (En) = e-[-^(^0,nA^)-r,(i?o,(n-l)A.)] ^ ^^^^^^ _ ^^29) 

Now, since Az ~ 10~^ evolutionary effects can be neglected inside a single domain and only 
accounted for when jumping from one domain to the next. As a consequence - owing to Eq. 
dTSl) - Eq. (^iM reduces to 



<l>(EO=e-^--/'^"'(^°)$(i?n+i) 



(130) 
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and the comparison of Eqs. (129) and (130) yields 



(n) 
dom 



{Eo, n Az) — Tj {Eq, {n — 1) A; 



(131) 



Further, by inserting Eq. (115) into Eq. (131), we get the desired photon mean free path 

Az 



4.29 • lO-' 



1 + 1.45 {n - l)Az J {Eo,nAz) - {Eq, {n - l)Az) 



Mpc . (132) 



F. Propagation over many domains 

We are finally ready to carry the strategy outlined in Subsection IV-B to completion, 
namely to evaluate the photon survival probability P^_^^^^{Eq, z) for a monochromatic 
beam emitted by a blazar at redshift z and detected at energy Eq. 

This task can be accomplished by first noticing that for a considered blazar at redshift z 
the overall behaviour of the photon/ALP beam is described by the following transfer matrix 



(133) 



n=l 



According to Eq. (104), the probability that a photon/ALP beam emitted by a blazar at z 



in the state pi will be detected in the state p2 for fixed orientations ipi, ...jipNa of B in every 
domain is 



P, 



Pl-S>P2 



{Eo,z;'^i,...,'^Na) = Tr (p2W(Eo,z;V^i,...,V^7vJ Pi W"^ (^o, ^; ^i, ^jvj) , (134) 



where it is assumed that Trpi = Trp2 = 1. As a consequence, the actual detection probability 
for the beam in question emerges by averaging the above expression over all angles, namely 



p 



{Eo,z;ipi,...,ipNa) 



(135) 



Because of the fact that the photon polarization cannot be measured at the energies consid- 
ered here we have to sum this result over the two final polarization states 



(136) 





p. = I 1 




(137) 



Moreover, we suppose for simplicity that the emitted beam consists 100 % of unpolarized 
photons, so that the initial beam state is described by 



Punpol 




(138) 
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M/(10ii GeV) 


5o/nG 


0.1 


1 - 60 


0.1 - 6 


0.5 


1 - 12 


0.5 ~ 6 


1.0 


1 - 6.0 


1 - 6 


5.0 


1 - 1.2 


5-6 



TABLE II: Allowed values of M and Bq in the considered cases. 



Hence, we ultimately have 

pDARMA (^^^^ 



7^7 V"U5~y \ " Punpol^'Pa: 



iJi,...,^N, 



+ 



(139) 



We implement this procedure as follows. In the first place, we arbitrarily choose the 
angle ipn in the n-th domain and we evaluate the corresponding transfer matrix W„ {En, ipn) 
for a given value of Eq, keeping Eq. (123) in mind. Next, the application of Eqs. ( |133 ) 
and (134) yields the corresponding photon survival probabilities entering Eq. ( 139[ ) for a 
single realization of the propagation process. We repeat these steps 5000 times, by ran- 
domly varying all angles ipn each time, thereby generating 5000 random realizations of the 
propagation process. Finally, we average the resulting photon survival probabilities over all 
these realizations of the propagation process, thereby accomplishing the average process in 
Eq. (139). We find in this way the physical photon survival probability P^^^^ {Eo,z). 



V. DISCUSSION 



Let us proceed to investigate the implications of the DARMA scenario for VHE blazars 
observations. 

We begin by stressing that all its physical predictions depend solely on B/M and not 
on B and M separately (this was true for photon-ALP oscillations and it remains true in 
general, because absorption does not depend on these quantities). For this reason, it is quite 
useful to introduce the dimensionless parameter 



So 
nG 



10" GeV\ 
M ) ' 



(140) 



Owing to conditions (60) and (111), it will be assumed 



e<6 



(141) 



throughout our discussion. Specifically, we will focus our attention on the representative 
cases ^ = 5.0, ^ = 1.0, ^ = 0.5, ^ = 0.1, taking both Ldom = 4Mpc and Ldom = lOMpc 
at z = 0. Nevertheless, it is important to keep under control which values of Bq and M 
are allowed in each case. From the constraints (60) and (111) we find the allowed ranges 
reported in Table [11} 

Next, we have to make sure that we stay within the strong-coupling regime all the way 
up to the source for Eq > 100 GeV. Therefore, by combining Eq. (44) with the requirement 
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Upper bound on m 


Value of the parameter 


4.40 • 10-1° eV 


C = 5.0 


1.97- 10-1° eV 


e = i.o 


1.39 • 10-1° eV 


C = 0.5 


0.62 • 10-1° eV 


C = 0.1 



TABLE III: Upper bounds on the ALP mass in the considered cases. 



< 100 GeV the resulting upper bound on m can be expressed as 



Im^ -cj2jV2 < 1.97. 10 



-10 



B 



1/2 



10iiGeV\'/' 



(142) 



which - on account of Eqs. (107) and (140) - can be more suitably rewritten in the form 

1/2 



m 



10-10 eV 



(1.14-10-^)' 



< 1.97^^/2 



(143) 



valid for all sources considered here |121] . Thanks to condition (141 ), we see that within the 
DARMA scenario ALPs have to be very light, with mass never exceeding 5 ■ 10-i°eV. In 
particular, the axion needed to solve the strong CP problem is therefore ruled out by several 
orders of magnitude. Observe that for m < 1.14 ■ IQ-i^eV the plasma frequency dominates, 
so that even massless ALPs behave as if their mass where equal to the plasma frequency. 
The upper bounds on m corresponding to the cases under consideration are reported in 
Table [Till 

As far as EBL absorption is concerned, we will take for the optical depth entering Eq. 

provided by the FRV model [97] . 
because in the absence of EBL absorption photon-ALP 
an enhancement of the photon survival probability 



(132) the exact expression {Eq, z) 



A general expectation is that 
oscillations only produce a dimming 
with respect to the case of conventional physics shows up only at sufficiently high energy, 
where EBL absorption becomes substantial. Therefore, close enough to 100 GeV a dimming 
rather than an enhancement should occur. 



VI. PREDICTIONS FOR FUTURE OBSERVATIONS 

The best way to figure out the relevance of the DARMA scenario for future observations 
to be performed with the CTA and with the HAWC water Cherenkov 7-ray observatory 
is to compare the photon survival probability P]^_^^^^{Eq, z) with the one predicted by 
conventional physics P^E,^{Eq, z), with the EBL described in either case by the FRV model. 

We do that for a sample of different redshifts, like z = 0.031, z = 0.188, z = 0.444 and 
z = 0.536. We remark that the case of z = 0.031 may look somewhat academic, since 
its location inside the Local Group is likely to make the morphology of the magnetic field 
crossed by its line of sight more complicated than assumed in this paper. Nevertheless, we 
include z = 0.031 in the present analysis in order to see what happens for a very nearby 
blazar even if a drastic simplifying assumption is made. 

The results are displayed in Figure 3. For each of the selected sources, we consider the 
above choices for ^, which are represented by a solid black line = 5.0), a dotted-dashed 
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line = 1.0), a dashed line = 0.5) and a dotted line = 0.1), while the solid grey line 
corresponds to conventional physics. We take both Ljom = 4Mpc and Ldom = lOMpc for 
the domain size at z = 0. 

All plots show one common trend. At energies only slightly in excess of 100 GeV, 
P^X-yi-^o, is larger than P^:^^^^{Eq, z), indeed in agreement with expectations. As the 
energy further increases, the situation reverses and P^_^^^^{Eq, z) gets progressively larger 
and larger than P^^(£'o, z) until the value 100 TeV is attained which is the highest energy 
value considered in the present analysis. 

A somewhat surprising result emerges at large enough energies. Indeed, since ^ sets 
the strength of the photon-ALP oscillation mechanism, it would be natural to expect 
P^-^^^^{,Eq, z) to monotonically increase with ^. However, this is not the case. More 
specifically, for ^ = 5.0 the behaviour of P^^^^^{Eq, z) as a function of Eq resembles 
closely that of P^E,^{Eq, z) - apart from an overall shift towards higher energies - and it 
is practically independent of Ldom apart from the case of 2; = 0.031 which exhibits a mild 
-Z^dom-dependence. Moreover, at sufficiently high energies the values of P^:^^^^{Eq, z) cor- 
responding to ^ = 5.0 are the lowest predicted by the DARMA scenario for all sources. 
The case ^ = 1.0 is different, since the resulting values of PJ^^^^^{Eq, z) are among the 
highest predicted at low redshift but become among the lowest at high redshift even though 
they always exceed those corresponding to ^ = 5.0. In addition - with the exception of 
z = 0.031 - PI^^J^^^{Eo, z) strongly increases as Ldom decreases. As ^ decreases the trend 
exhibits a radical modification. Actually, the case ^ = 0.5 shows a mild Ldom-dependence 
for all considered sources, and with the exception of z = 0.031 it leads to the largest values 
of P^_^^^^{Eo, z) for Ldom = 4Mpc. Finally, the case ^ = 0.1 depends more strongly on 
Ldom and - again with the exception of z = 0.031 - for some energy values it can make 
P^^^^^{Eq, z) larger than in the case ^ = 0.5 for Ldom = 10 Mpc, but this never occurs for 
Ldom = 4 Mpc. We observe that the different situation found for z = 0.031 as compared to 
the other blazars should not come as a surprise, owing to the above remarks. 

What is the reason for such a behaviour? 

Owing to the random structure of the considered magnetic field, coherence is maintained 
only within one domain and so P}^_^^^^{Eo, z) is ultimately controlled by two quantities: 
the photon-ALP conversion probability over a single domain P^^a{Ldom) and the photon 



absorption probability in Eq. (84). In order to clarify this issue in an intuitive fashion, we 



argue as follows, discarding cosmological effects for simplicity. 

As far as P-y^ai^L^o^ is concerned, we have seen that it is given by Eqs. (42) in the case 
of photons linearly polarized in the direction parallel to B. This is not true in the present 
situation where the beam photons are assumed to be unpolarized, but for the sake of an 
order-of-magnitude estimate we can still suppose that -Ry-s.a(Ldom) has the form (42) and 
therefore we write it as 



P^^a (Ldom) - sin^ 



1.6-10""^ 



Mpc 

We distinguish two cases and we discuss them in turn: 



(144) 



As long as ^ <^ 60 (Mpc/Ldom), Eq- (144) yields Py->,a(Ldom) ^ 1 which entails 



that the fraction of ALPs produced over a single domain is very small. Since we 
are supposing the beam to be initially fully made of photons, it takes a length much 
larger that Ldom before a sizeable fraction of the beam consists of ALPs. In the same 
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FIG. 3: Behaviour of PJ^^^^^ versus the observed energy Eq for: z = 0.031 (top row), z = 0.188 
(second row), z = 0.444 (third row), z = 0.536 (bottom row). The sohd black Hne corresponds to 
^ = 5.0, the dotted-dashed Hne to ^ = 1.0, the dashed hne to ^ = 0.5, the dotted hne to ^ = 0.1 
and the sohd grey hne to conventional physics. We have taken Ldom = 4Mpc (left column) and 
-^^dom = lOMpc (right column). 
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FIG. 4: Behaviour of P^^^^^ versus the observed energy Eq for: z = 0.031 (top row left), 
z = 0.188 (top row right), z = 0.444 (bottom row left), z = 0.536 (bottom row right) for Ldom = 
0.05 Mpc. The solid black line corresponds to ^ = 5.0, the dotted-dashed line to ^ = 1.0, the 
dashed line to ^ = 0.5, the dotted line to ^ = 0.1 and the solid grey line to conventional physics. 



fashion, once such a situation is realized, a similar long length is needed in order for 



the beam to contain a sizeable amount of photons. Moreover, it follows from Eq. (144) 
that presently -P^^-a (-^dom) becomes a quadratic function of ^ which therefore increases 
monotonically with ^. Accordingly, the picture outlined in Subsection IV-A is expected 
to emerge straightforwardly and this is confirmed by a numerical simulation in which 
^ takes the above values but we assume Ldom = 0.05 Mpc, which yields the behaviour 
shown in the plots reported in Figure 4. 

When condition ,^ ^ 60 (Mpc/Ldom) is not fulfilled the situation becomes considerably 
more complicated. In the first place, P^-^a{Ldom) fails to be a monotonically increasing 
function of C, and it becomes oscillatory. So, depending on the actual value of Ldom 



it follows that P^^aiLdom) can decrease as ^ increases. As stressed above, Eq. (144) 
can be taken at most to provide an order-of-magnitude estimate but it is clear that 
condition ,^ ^ 60 (Mpc/Ldom) fails to be met for L^om = 4 Mpc and L^om = 10 Mpc 
along with the considered values of ^. Hence, for a fixed source distance the behaviour 
exhibited in the plots in Figure 3 can arise. Still, this is not the end of the story, since 
in the present situation even after the domain closest to the source a relevant fraction 
of the beam consists of ALPs. In other words, a large enough number of 7 — )■ a and 
a — )■ 7 transitions take place inside a single domain. So, the overall effect is to have a 
larger number of photons per unit length between the source and us as compared to the 
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previous case. As a consequence, EBL absorption is now more effective thereby giving 
rise to a smaller observed photon ffux. Moreover, this dimming evidently increases 
with the source distance, which explains why P^^J^^^{Eq, z) tends to decrease as z 
increases at fixed energy and eventually it behaves like P^]l^{Eo, z), which indeed 
occurs for ^ = 5.0. 

In conclusion, it is evident from Figures 3 that in the most favourable case a boost 
factor of 10 in P}^_^^^^{Eo, z) occurs at progressively lower energies as the source distance 
increases. Explicitly, for z = 0.031, z = 0.188, z = 0.444, z = 0.536 it takes place at 
EiQ ~ 30 TeV, 8 TeV, 2 TeV, 1.2 TeV, respectively. Above Eio the boost factor can be much 
larger. 



VII. A NEW INTERPRETATION OF OBSERVED VHE BLAZARS 

Observed VHE blazars provide a great deal of information which can challenge the 
DARMA scenario. Within conventional physics the values of Fem have to be tuned for 
every source in such a way to reproduce the corresponding values of Fobs- No rational lies 
behind this procedure and a large spread in the values of Fem is demanded in order to account 
for the equally large spread in the values of Fobs- While this procedure does not pose any 
technical problem, such a systematic fine-tuning lacks any conceptual appeal and moreover 
leads to the cosmic opacity problem. 

Thus, it looks natural to inquire whether the DARMA scenario sheds some light on this 
issue. 

The most straightforward way to investigate this question is to proceed somehow in 
parallel with the treatment followed in Subsection III-C, namely to de-absorb the values of 
Fobs within the present context. 

Our starting point is the general relation between the observed and emitted fluxes ex- 



pressed by Eq. (64), which presently reads 



$obs(i?0, Z) = P°™^(i?0, $ em {Eo{l + z)) . (145) 



Thanks to Eq. (77), we first rewrite Eq. (145) as 



K Ef"-^^'^ = P°™^(i?o, Z) $em {Eo{l + z)) , (146) 

from which we get $em(-E'o(l + z)) for every detected VHE blazar. We next best-fit this 



function to the power-law expression (76), namely 



_pDARMA 

^,^{Eo{1 + z)) = k\Eo{1 + z)\ (147) 

over the energy range where the considered source is observed (see Table j!)). We obtain in 
this way the values of F^j^^"^^. 

We stress that for a given choice of ^ and Ldom the photon survival probability 
-P^^'^^(-E'o, z) is uniquely fixed, apart from errors affecting r^^^(£'o, z) which are unknown 
and therefore again ignored. Since F^^^^ is linearly related to Fobs? the associated error 
bars are the same even in the present context (see Table [l]) and so they will not be explicitly 
exhibited. 
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We implement the considered procedure by taking for the free parameters ^ and Ldom 
the same representative values chosen above. The resulting numerical values of F^^^"^ for 
the various cases are reported in Tables IV and |V| together with the corresponding spread 
AF^^^"^ of the value of Fcm neglecting errors and the average value (F^^'^'^) including 
errors. 

Let us consider first what happens as ^ increases from ^ = 0.1 to ^ = 5.0 assuming 
-^^dom = 4Mpc. As long as ^ = 0.1 the difference between the DARMA scenario and 
conventional physics strictly vanishes up to 2; = 0.138 and becomes negligible at larger 
redshifts. However, as soon as the regime C, = 0.5 is attained the DARMA scenario starts 
to differ drastically from conventional physics. As ^ increases from 0.5 to 5.0 it is found 
that (F^^^^) monotonically increases, even if at a rate that slows down for increasing ^, 
which entails that no improvement is to be expected for ^ > 5.0 regardless of any other 
consideration. On the other hand, AF^^'^^ decreases for ^ = 0.1— 7-^ = 1. but next 
increases for ^ = 1.0 — )■ ^ = 5.0. 

A somewhat similar pattern is found for L^om = 10 Mpc. Again for ^ = 0.1 the difference 
between the DARMA scenario and conventional physics strictly vanishes up to 2; = 0.116 
and remains negligible at larger redshifts. A big difference shows up around ^ = 0.5. As 
before, (F^^^^^^) increases monotonically, but its rate slows down for increasing ^ and just 
vanishes in the step ^ = 1.0 — )■ ^ = 5.0. Here AF^^^"^ still decreases for ^ = 0.1 — )■ ^ = 0.5 
but then increases for^ = 0.5—)'^ = 5.0. 

Let us next find out what happens in the change Ldom = 4 Mpc — L^om = 10 Mpc at 
fixed ^. As far as AF^^'^^ is concerned, it decreses for ^ = 0.1 and ^ = 0.5 but it increases 
for ^ = 1.0 while it remains practically unchanged for ^ = 5.0. The behaviour of (F^^^^^^) 
is slightly different, since only in the case ^ = 0.5 it shows a slight variation. 

Physically, all this means that for ^ = 0.1 DARMA effects are negligible, but they sud- 
denly become important shortly before ^ = 0.5 is reached and they remain more or less 
unchanged up to ^ = 5.0. This conclusion is in remarkable agreement with our previous 
results concerning the behaviour of P}^_^^^^{Eq, z). In this connection, two points should 
be stressed. It follows from Eqs. (146) and (147) that at fixed z F^^^^ depends logarith- 
mically on P}^_^^^^{Eo, z), which makes its dependence on ^, Eq and Ldom niuch shallower 
than that of P^^^^(-Eo, z) itself. Consequently, the sharp differences found in Section VI 
among the cases with different ^ get smoothed out here, apart from one thing: since in 
the case ^ = 0.1 the photon survival probability is nearly the same in conventional physics 
and within the DARMA scenario over the energy range 0.2 TeV < Eq < 2TeV where most 
blazars are observed, the same is evidently true for the values of Fem- In addition, the 
marked difference among the cases ^ = 0.1, 0.5, 1.0, 5.0 discovered in Section VI takes place 
at energies considerably larger than those in the presently considered range, which explains 
why the cases ^ = 0.5, 1.0, 5.0 exhibit a fairly similar behaviour for the observed blazars. 



A. Solution of the cosmic opacity problem 

A glance at Tables IV and [v] shows that the values of F^^'^^ for all VHE blazars 
happen to be in the same ballpark, thereby implying that within the DARMA scenario the 
observations can be explained with the same physical mechanism operating in all blazars - 
as a consequence, there is no cosmic opacity problem. 
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Clnnrpp r^ARMA pDARMA -pDARMA -p 

ouuiue ■■- em ^ cm ^ cm -"-em 



^ = 0.1 ^ = 0.5 ^ = 1.0 ^ = 5.0 



3C 66B 


3.00 


3.00 


3.00 


3.00 


3.03 


Mrk 421 


2.16 


2.16 


2.16 


2.17 


2.21 


Mrk 501 


1.90 


1.90 


1.90 


1.91 


1.96 


Mrk 501 


2.03 


2.03 


2.03 


2.04 


2.08 


lES 2344+514 


2.70 


2.70 


2.71 


2.73 


2.78 


Mrk 180 


3.07 


3.07 


3.07 


3.09 


3.14 


lES 1959+650 


2.43 


2.43 


2.44 


2.46 


2.53 


BL Lacertae 


3.27 


3.27 


3.28 


3.32 


3.38 


PKS 0548-322 


2.39 


2.39 


2.40 


2.45 


2.52 


PKS 2005-489 


3.59 


3.59 


3.60 


3.66 


3.73 


RGB J0152+017 


2.47 


2.47 


2.48 


2.56 


2.63 


W Comae 


3.18 


3.18 


3.21 


3.32 


3.39 


PKS 2155-304 


2.67 


2.67 


2.72 


2.85 


2.90 


RGB J0710+591 


? 


? 


? 


? 


? 


H 1426+428 


0.85 


0.85 


1.28 


1.57 


1.44 


lES 0806+524 


2.70 


2.70 


2.77 


2.93 


3.00 


lES 0229+200 


0.41 


0.42 


1.15 


1.37 


1.13 


H 2356-309 


2.06 


2.06 


2.17 


2.35 


2.40 


lES 1218+304 


2.00 


2.00 


2.15 


2.32 


2.36 


lES 1101-232 


1.72 


1.73 


1.96 


2.13 


2.13 


lES 0347-121 


1.87 


1.87 


2.11 


2.28 


2.28 


lES 1011+496 


2.90 


2.90 


3.06 


3.22 


3.26 


S5 0716+714 


1.60 


1.61 


2.07 


2.22 


2.22 


PG 1553+113 


2.48 


2.49 


3.00 


3.08 


3.08 


PKS 1222+21 


2.47 


2.47 


2.80 


2.88 


2.90 


3C 66A 


1.28 


1.30 


2.19 


2.25 


2.22 


PKS 1424+240 


1.16 


1.18 


2.03 


2.06 


2.04 


3C 279 


2.05 


2.06 


2.71 


2.74 


2.73 




3.18 


3.17 


2.45 


2.29 


2.60 


(Pem) 


2.22 


2.23 


2.41 


2.51 


2.52 



TABLE IV: We have inferred the emitted spectral index Pem by de-absorbing within the DARMA 
scenario the observed value of Pobs for every source neglecting errors. This proccdTirc has been 
carried out for the choice of parameters ^ = 0.1, ^ = 0.5, ^ = 1.0 and ^ = 5.0. In all cases, we 
have taken L^om = 4Mpc. The similar values obtained in Subsection III-C within conventional 
physics have been quoted for comparison and are denoted by T^^. The last two lines report the 
spread APem of the value of Pem discarding errors and the average value (Pem) including errors, 
respectively, for the various cases. 

B. Fitting individual sources 

The foregoing analysis has shown that the gist of the DARMA scenario for the observed 
VHE blazars is to drastically reduce the spread in the values of Pem as compared with what 
happens in conventional physics, thereby tracing the large spread in the values of Fobs to 
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Source 



pCP pDARMA 

em em 



■pDARMA 

em 



■pDARMA 

em 



■-DARMA 









e = 


= 0.1 


e = 


= 0.5 


e = 


= 1.0 


^ = 5.0 


3C 66B 


3 


00 


3 


00 


3 


03 


3 


01 


3.03 


Mrk 421 


2 


16 


2 


16 


2 


16 


2 


18 


2.21 


Mrk 501 


1 


90 


1 


90 


1 


90 


1 


93 


1.96 


Mrk 501 


2 


03 


2 


03 


2 


03 


2 


06 


2.08 


lES 2344+514 


2 


70 


2 


70 


2 


72 


2 


76 


2.79 


Mrk 180 


3 


07 


3 


07 


3 


08 


3 


11 


3.14 


lES 1959+650 


2 


43 


2 


43 


2 


45 


2 


50 


2.53 


BL Lacertae 


3 


27 


3 


27 


3 


30 


3 


36 


3.38 


PKS 0548-322 


2 


39 


2 


39 


2 


43 


2 


50 


2.52 


PKS 2005-489 


3 


59 


3 


59 


3 


63 


3 


70 


3.73 


RGB J0152+017 


2 


47 


2 


47 


2 


53 


2 


61 


2.63 


W Comae 


3 


18 


3 


18 


3 


28 


3 


37 


3.39 


PKS 2155-304 


2 


67 


2 


67 


2 


81 


2 


89 


2.90 


RGB J0710+591 




? 




? 




? 




? 


? 


H 1426+428 





85 





86 


1 


53 


1 


53 


1.44 


lES 0806+524 


2 


70 


2 


71 


2 


88 


2 


98 


3.00 


lES 0229+200 
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48 


1 


36 


1 


27 


1.12 


H 2356-309 


2 


06 


2 


06 


2 


31 


2 


39 


2.40 


lES 1218+304 


2 


00 


2 


00 


2 


29 


2 


35 


2.36 


lES 1101-232 


1 


72 


1 


73 


2 


10 


2 


14 


2.13 


lES 0347-121 


1 


87 


1 


88 


2 


26 


2 


29 


2.28 


lES 1011+496 


2 


90 


2 


90 


3 


19 


3 


25 


3.26 


S5 0716+714 


1 


60 


1 


62 


2 


20 


2 


22 


2.22 


PG 1553+113 


2 


48 


2 


52 


3 


07 


3 


08 


3.08 


PKS 1222+21 


2 


47 


2 


49 


2 


87 


2 


89 


2.90 


3C 66A 


1 


28 


1 


36 


2 


25 


2 


23 


2.22 


PKS 1424+240 


1 


16 


1 


26 


2 


07 


2 


05 


2.04 


3C 279 


2 


05 


2 


13 


2 


74 


2 


74 


2.73 




3 


18 


3 


11 


2 


27 


2 


43 


2.61 


(Pem) 


2 


22 


2 


24 


2 


49 


2 


52 


2.52 



TABLE V: Same as Table IV but with Ldom = 10 Mpc 



the wide spread in the blazar distances. 

It seems therefore worthwhile to investigate this point in a quantitative fashion according 
the following strategy: 

• As a zero-order approximation, we suppose that all blazars have the same value of 
r^^'^^, which for definiteness is taken to be the average value over all observed 
sources (F^^^^) for a given choice of ^ and Ldom- 

• As a first-order correction - which is meant to improve on the above idealized situation 
- we allow for a small spread around (F^^^^), which we tentatively take to be ± 0.2. 

In order to keep the situation under control, we focus our attention on the single case 
^ = 1.0 and Ldom = 4 Mpc which we regard as the most favourable one not only because 
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Ar^^'^^ is very small - the case ^ = 0.5 and Ldom = 10 Mpc would be even better in this 
respect - but also because we feel that L^om = 4 Mpc is more realistic than Ldom = 10 Mpc. 
This amounts to take (r°^^^^) = 2.51, which entails in turn 2.31 < r°^^^^ < 2.71 for all 
observed VHE blazars. The value 2.51 is close to the value 2.40 that we used in a previous 
discussion of the DARMA scenario [25J, as well as to 2.47 which is the average value for the 
observed VHE blazars with z < 0.05 that undergo a negligible EBL attenuation. 

Next, we evaluate for every source the e3;]:>ecie(^ observed spectral index T'^^{z). Basically, 
this amounts to run backwards the same procedure whereby we have got the values of 



■^DARMA 



reported in Table IV for = 1.0. Explicitly, by combining Eqs. (145) and (147) we 



can write 



pDARMA/ 77. 



K 



Eo{l + z) 



-2.51 



(148) 



Since P^^^^{Eq,z) is known, ^^^^(i^o,-^) can be computed exactly. Then we best-fit this 



function to the power-law expression (77), namely 



(149) 



over the energy range where each source is observed. We find in this way the values of 
-'^obs(^) every source. As repeatedly stressed, the observed and emitted spectral indices 
are linearly related, and so they have the same error bars. 

We are now ready to check this view by performing a fit to all observed VHE blazars. 
This is shown in Figures^ to ^ where the solid black line corresponds to r^^g(z) while the 



grey strip represents the range ^''^J^z) ± 0.2. We stress that r|^^g(z) is different for different 
sources, owing to the different observed energy range. 

A look at those Figures shows that by assuming that all VHE blazars have Fem in a 
range of 2.51 ± 0.2 allows to fit observations of 19 sources out of a total of 27 ones. The 
role of photon-ALP oscillations is to partially offset EBL absorption, and thus the DARMA 
scenario departs from conventional physics only to the extent that EBL attenuation becomes 
important; for 2; > 0.1 12 sources out of a total of 16 ones are successfully fitted, and for 
z > 0.138 the fit is successful for 10 blazars out of a total of 11. 



VIII. CONCLUSIONS 

Very light ALPs are a generic prediction of many attempts to extend the Standard Model 
along different directions towards a more satisfactory fundamental theory of all elementary- 
particle interactions including gravity. We have systematically investigated the DARMA 
scenario in which the mechanism of photon-ALP oscillations triggered by large-scale mag- 
netic fields is regarded as a means to effectively reduce the EBL attenuation affecting blazar 
observations above 100 GeV. Our assumptions can be summarized as follows: 

• Large-scale magnetic fields exist with a cellular morphology characterized by a coher- 
ence lenth in the 1 — 10 Mpc range and a strength not much smaller than the available 
upper bound i?o < 6 nG. 

• ALPs have to be very light in order to ensure that the strong-mixing regime is realized. 
The upper bound on their mass depends on the adopted value of the 077 coupling 



constant Bq/M - see Table HI - but in any case the condition m < 5 ■ 10~^° eV has to 
be met. This prevents the axion needed to solve the strong CP problem from playing 
any role in the present context. 
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• The parameter M is consistent but fairly close to the strongest upper bound M > 
10^^ GeV coming from observations of supernova SN1987a. We remark that this bound 
is however affected by a large uncertainty and exceeds by one order of magnitude the 
robust bound M > 10^° GeV coming both from theoretical considerations of star 
cooling and from the negative result of the CAST collaboration. 

We predict that a boost factor of 10 in the photon survival probability with respect to 
conventional physics takes place for all VHE blazars observed so far well below the upper 
detection threshold of the planned CTA and HAWC water Cherenkov 7-ray observatory. 
Moreover, the energy Eiq at which such a boost factor occurs decreases as the source distance 
increases and becomes e.g. as low as 2TeV for the blazar 3C 279 at z = 0.536. Hence, our 
prediction can certainly be tested with the above planned detectors and possibly also with 
currently operating lACTs H.E.S.S., MAGIC, CANGAROOIII, VERITAS as well as with 
the Extensive Air Shower arrays ARGO-YBJ and MILAGRO. 

We find it a remarkable fact that the DARMA scenario also offers a new interpretation of 
the observed VHE blazars, according to which the values of Fem for far-away VHE blazars 
are in the same ballpark of nearby ones and the large spread in the values of Fobs is mainly 
traced to the wide spread in the source distances. 

As is well known, weakly interacting massive particles (WIMPs) can be detected either 
indirectly through astrophysical effects or directly at the Large Hadron Collider (LHC). The 
situation of ALPs characteristic of the DARMA scenario is in a sense similar. Besides being 
detectable indirectly through the astrophysical effects discussed in this paper, they lend 
themselves to a direct detection either in the GammeV [l3j experiment at FERMILAB or 
more likely with planned photon regeneration experiment ALPS at DESY [H], or else with 
large xenon scintillation detectors developed for dark matter searches |45j . 

Note added in proof. After submission of the present paper, we have become aware 
that our predicted lower-than-expected transparency of the Universe in the VHE band is 
supported by an independent result |124] which rests upon a new statistical analysis of all 
VHE blazars based on a Kolmogorov-Smirnov test in conjunction with the minimal EBL 
model 1921 . 



Many enlightening conversations with our collaborators Oriana Mansutti and Massimo 
Persic are gratefully acknowledged. We also thank Dieter Horns and Fabrizio Tavecchio for 
informative discussions and Aldo Treves for several remarks. 



We solve here the mathematical problem of finding the transfer matrix U{y, i/q] 0) asso- 
ciated with the reduced Schodinger-like equation 
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(150) 
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(151) 




(152) 



X2 = \ [ t + u - J {t - uf + Av"^ 



with 

/ My) 

as in the text, and mixing matrix of the form 



M = 

where the coefficients s, t, u and v are supposed to be complex numbers. 
We start by diagonahzing M.. Its eigenvalues are 

Ai = s , (153) 

(154) 

(155) 

and it is straightforward to check that the corresponding eigenvectors can be taken to be 

Xi = I I , (156) 

(157) 
(158) 

(159) 




Correspondingly, any solution of Eq. (150) can be represented in the form 



i){y) = ci Xi e'^' + C2 X2 e'^^ (^"^"^ + C3 X3 e'^' ^^-^'^ 



where Ci, C2, C3 and yo are arbitrary constants. As a consequence, the solution with initial 
condition 

Ax{yo) 

ij{yo) = I AM) I (160) 

a{yo) 



emerges from Eq. (159) for 



C2 



ci = A^{yo) 
Xs-t . , . 1 



C3 = - 



'^(Aa - A2) 
A2-t 



Aiyo 



v{X3 - A2) 



A{yo) + 



A3 — A2 
1 



A3 -A 



a{yo) , 
-a{yo) ■ 



(161) 
(162) 

(163) 
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It is a simple exercize to recast the considered solution into the form 

^(l/)=W(i/,yo;0)V^(yo) (164) 

with 

U(y^ y^- 0) = e*^i(f-f") Ti(0) + e'^^^y-y°^ T2(0) + e^^«(^-^«) r3(0) , (165) 
where we have set 

/ 1 0\ 

(166) 



T2(0)^ , aS , , (167) 



/I 0^ 




= 000 




\o oy 










\i-t 


V 


A3— A2 

(A2-t){A3-t) 


A3— A2 

A2-t 


'U{A3-A2) 


A3-A2 








A2-t 


V 


A3— A2 

(A2-t)(A3-t) 


A3— A2 


'U{A3-A2) 


A3-A2 



T3(0)^ , aA , (168) 



from which it follows that the desired transfer matrix is just U{y,yo;0) as given by Eq. 



(165). 



X. APPENDIX B 

It proves very useful for illustrative purposes to have the approximate behaviour of 
T.y{E(),z) in an analytic form. This goal can be achieved by taking advantage from the 



fact that aj^{E,e,(p) is maximized when condition (70) is met (we restrict ourselves to 
head-on collisions for simplicity). 

Accordingly, the crudest attempt to estimate the dominant contribution to the optical 
depth would be to approximate the e{z) integration by the product of a^^"^ times n^{e{z),z) 



as evaluated for that particular value of e{z) selected by condition (70) for fixed E. This 
amounts to insert the Dirac delta 6{e{z')/eV — 500 GeV/i?(2;')) into the r.h.s. of Eq. (71), 
which leads to 



T^{Eo, Z) ~ 2.25 • 10^ / AZ '"^ ^^i ^ (^-^gg^i 



^7 (e(^),^) 







l + z) [0.7 + 0.3(1 + 2)^] 



31 1/2 



with 

1 /500GeV\ 
e(.). — jeV. (170) 



Unfortunately, experience with this problem shows that the resulting £'o-dependence of 
T^{Eo, z) is too steep, and since P^^(i?o, z) depends exponentially on t^{Eq, z) this approx- 
imation is doomed to failure. 

A more satisfactory conclusion emerges by exploiting a popular approximation [122] 
which amounts to replace the e{z) integration in Eq. (71) by the product of a'^J^^ times 
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e{z) ny{e{z), z) at the particular value of e{z) dictated by condition (70) for fixed E. We 
find in this way 



T^{Eo,z) ~ 2.25 • 10^ 



500GeV\ 



dz 



[l + zY [0.7 + 0.3(1 + ^)'] 



31 1/2 



(171) 



(500 GeV/Ep) eV 



z cm eV 



where Eq. (170) has been used 



Owing to Eq. (171), the derivation of an approximate analytic behaviour of the optical 
depth requires an approximate analytic expression for n^{e{z),z). Unfortunately, the FRV 
model does not give an analytic form for n^{eo, 0) but provides a plot of eo n^(eo, 0) versus eo 
(see their Figure 4). The corresponding plot of n^(eo, 0) as a function of eo is reproduced by 
the solid line in Figure [TO] for the EBL energy range relevant for the observed blazas, namely 



0.25 eV < Co < 2.5 eV. A look at Figure 10 shows that within the considered energy range 
the spectral energy distribution of the EBL departs from a power-law behaviour, owing to 
the emission bump resulting from the integrated emission of the low-mass stellar population 
that remained close to the main-sequence over cosmological times. For this reason, we 
approximate the FRV result for n^(eo, 0) with the shadowed linear strip shown in Figure 10 
enveloping the exact behaviour, which is expressed by the following power-law representation 



n] 



app. 



eo,0) 



10"^ a 



eV\ 



1.85 



cm eV 



(172) 



with the constant a in the range 0.9 < a < 3.6 so as to enclose the FRV curve. 

Evolutionary effects in the EBL spectral energy distribution can be taken into account as 
follows. Besides redshifting all energies in proportion of 1 -|- 2;, the cosmic expansion dilutes 
the EBL by a factor (1-1-2;)^ similarly to what happens for the CMB. But in addition the EBL 
spectral photon number density changes because of the intrinsic evolution of the galactic 
population over cosmic times. A quantitative analysis [123J shows that the EBL photon 
number density acquires an extra factor (1 -|- z)'^'"^ as long as z < 1, which is appropriate 
to our case. On the whole, the spectral photon number density n'^^^{e{z), z) of the EBL at 
redshift z is related to n^^^{eo, 0) by 



<PP(e(z), z) deiz) ~ (1 + z)'-' <PP(eo, 0) deo 



which yields 
namely 



<PP(e(^),z)^(l + z)°%r(^,0 



n: 



app 



(500 GeV/Ep) eV 
l + ~z ' 



[l + z 



,0.8 ^app 



(500 GeV/Ep) eV 



(173) 

(174) 

(175) 



thanks to Eq. (|170|). In particular, Eq. (|172|) leads to 



app 



(500 GeV/Ep) eV 
l + ~z ' 



500 GeV 



1.85 



4.5 -3 -ir-l 

I cm eV 



(176) 
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The approximate evaluation of the optical depth t^^^^Eq, z) amounts to insert Eq. (176) 



into Eq. (171). Correspondingly, we get 

0.85 

500 GeV 



P(Eo, z) ^ 2.25 a ( — ^ ) I{z) , (177) 



where we have set 



I{z) ^ I dz' -— . (178) 

[0.7 + 0.3(1 + ^0 ] 



This integral has been evaluated numerically and its behaviour is depicted in Figure 11 
which shows that up to z ~ 0.1 it goes linearly with z but it increases more rapidly for 
larger redshifts. 



Before leaving this issue it is worthwhile to show how Eq. ( 73 ) is recovered in the limit of 
small z where cosmological effects become irrelevant. In such a situation the source distance 
is D = cz/Hq, which allows us to write the optical depth as 



E V' "" I D 



rr(£,«) = 5.25.10-a( ^l I ^ ) , (179) 



with the replacement i?o ^ So, we see that t^^^{E, D) oc D in agreement with Eq. (73), 
which entails that in the present approximation the mean free path for 77 — )■ e^e~ is given 

by 

A^PP(E) = 1.90 ■ 10=^ a-' ^ 500GeV \ ^^^^^ 



This quantity is plotted in Figure [12] as a function of E, where it is represented by the 
shadowed region between the dotted line corresponding to a = 0.9 and the dashed line 
corresponding to a = 3.6. The solid curve yields \-y{E) as evaluated exactly within the FRV 
model and exhibited in Figure [1} We see that the present approximation is indeed consistent 
with the result of the FRV model for 0.2 TeV <E <2 TeV. 
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FIG. 5: Behaviour of T^^^^"^ for the blazars 3C 66B, Mrk 421, Mrk 501 (with the two measure- 
ments of Fobs in the hterature), lES 2344+514 and Mrk 180. The sohd black hne corresponds to 
pDARMA ^ 2.51 and the grey strip represents the range 2.31 < T^j^^'^^ < 2.71. 
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FIG. 6: Behaviour of T^^™^ for the blazars lES 1959+650, BL Lacertae, PKS 0548-322, PKS 
2005-489, RGB J0152+017 and W Comae. The sohd black hne corresponds to r^^^^^^ = 2.51 
and the grey strip represents the range 2.31 < F^^'^^ < 2.71. 
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FIG. 7: Behaviour of T^^'^^^'^ for the blazars PKS 2155-304, H 1426+428, lES 0806+524, lES 
0229+200, H 2356-309 and lES 1218+304. The sohd black hne corresponds to r°j^^^^ = 2.51 
and the grey strip represents the range 2.31 < F^^'^^ < 2.71. 
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FIG. 8: Behaviour of r°^™^ for the blazars lES 1101-232, lES 0347-121, lES 1011+496, S5 
0716+714, PG 1553+113 and PKS 1222+21. The sohd black hne corresponds to T^^^^^^ = 2.51 
and the grey strip represents the range 2.31 < F^^'^'^ < 2.71. 




FIG. 9: Behaviour of TJ^^^^^'^^ for the blazars 3C 66A, PKS 1424+240 and 3C 279. The sohd black 
hne corresponds to F^^^^ = 2.51 and the grey strip represents the range 2.31 < T^^^^ < 2.71. 
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FIG. 10: The spectral photon number density in the present Universe n^(-,0) is plotted versus 
the energy cq in the energy range 0.25 — 2.5 eV. The sohd hne represents the result of the FRV 
model. The dotted and dashed lines correspond to the lower (a = 0.9) and upper (a = 3.6) limit, 
respectively, of our power- law approximation defined in Eq. (172). 
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FIG. 11: Plot of the behaviour of I{z). 
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FIG. 12: The approximate pair-production mean free path A^^^^ of a VHE photon is plotted versus 
its energy E and it is represented by the shadowed area as the parameter a varies in the range 0.9 
- 3.6. The dotted and dashed hues correspond to a = 0.9 and a = 3.6, respectively. Superimposed 
is the exact result obtained within the FRV model and shown in Figure [!} 



